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1.0  INTRODUCTION 


ModorJi  image  processing  systems  use  pattern  recognition  techniques 
to  automatically  process  imagery  generated  by  Forward  Looking  Infra-Red 
(FLIR)  and  other  sensors.  The  ATAC  FLIR  Autocuer  is  an  example  of  such 
a  system  that  employs  an  automatic  target  classifier.  The  initial  stage  in 
autocuers  is  the  extraction  of  shapes  from  the  scene  called  (object)  segmen¬ 
tation.  These  shapes  are  measured  and  the  resulting  features,  characteris¬ 
tic  of  targets,  are  input  to  a  classifier  for  automatic  target  classification/ 
recognition.  A  variety  of  algorithms  have  been  developed  to  perform  this 
shape  extraction  function.  These  algorithms  usually  employ  edge  detection 
and  intensity  thresholding  techniques,  and  some  type  of  connectivity  criteria 
which  are  heuristic  in  nature. 

The  performance  of  these  algorithms  has  been  judged  subjectively, 
usually  on  limited  data  bases  of  FLIR  imagery.  For  example,  segmentation 
algorithms  liave  been  termed  "good"  if,  when  applied  to  a  few  images,  the 
extracted  objects  look  good  to  the  algorithm  developers.  Factors  such  as 
noise  sensitivity,  presence  of  clutter,  or  pixel  quantization  are  usually  not 
a  part  of  the  evaluation  of  these  techniques,  even  though  these  effects  can 
have  a  profound  effect  on  the  performance  of  these  processes.  Quite  often, 
the  result  of  this  approach  is  a  set  of  algorithms  or  a  hardware  system  that 
performs  well  on  one  set  of  FLIR  images,  but  badly  on  another  set.  Then 
there  is  no  way  to  ascertain  qualitatively  the  reason  for  the  poor  perform¬ 
ance  or  just  what  should  be  done  to  improve  performance.  For  example,  it 
is  not  always  apparent  that  the  problems  are:  fundamental  to  the  algorithms; 
due  to  noisy  or  cluttered  data;  failure  of  the  hardware  or  software  imple¬ 
mentation;  or  just  a  statistical  accident. 


Useful  quauititative  (analytic)  results  relating  to  these  algorithms 
are  difficult  to  obtain,  due  to  the  difficulty  in  modeling  realistic  scenes,  the 
intractability  of  the  algorithms  themselves,  and  the  difficulty  in  defining 
relevant  and  meaningful  performance  parameters.  This  report  presents 
the  beginnings  of  a  more  formal  analysis  of  image  processing  systems  such 
as  the  ATAC  autocuer  and  in  particular,  the  segmentation  process.  Since 
the  autocuer  is  basically  designed  to  find  targets  in  a  FLIR  scene,  the 
application  of  these  results  must  be  considered  in  relation  to  the  following 
definitions. 


1. 

Detection 

Target  or  clutter 

2. 

Classification 

Class  of  target,  e.g.  , 
wheeled. 

tracked  or 

3. 

Recognition 

Type  of  vehicle,  e.g.. 

tank,  APC,  truck. 

4. 

Identification 

Type  of  tank 

The  goal  of  this  report  is  the  analysis  of  the  classification  of  targets,  with 
some  application  to  recognition.  In  particular,  many  of  the  current  results 
require  a  minimum  of  3  pixels  per  dimension.  The  results  will  apply  to 
recognition  to  ♦^he  extent  that  recognition  does  not  use  internal  structural 
detail,  that  is,  recognition  is  performed  by  shape  analysis  alone. 

The  models  generated  also  allow  for  simple  clutter,  where  clutter 
is  meant  to  be  anything  in  the  scene  that  is  not  a  target  or  part  of  a  uniform 
background.  Two  types  of  clutter  were  considered  in  terms  oi  their  effect 
on  segmentation  uf  a  nearby  target.  The  first  type  can  be  termed  object 
clutter,  such  as  a  rock  or  tree  that  exhibits  definite  shape.  The  second  type 
of  clutter  considered  could  be  termed  texture  in  that  it  does  not  exhibit 
a  definite  shape  but  rather  a  variation  in  the  background  intensity  in  the 
vicinity  of  a  target.  An  example  of  this  type  of  clutter  is  the  variation  in 
emissivity  due  to  soil  conditions. 

The  principal  results  of  this  analysis  can  be  summarized  as  follows. 
First  there  is  now  an  objective  qualitative  means  of  evaluating  shape  extrac¬ 
tion  algorithms  in  terms  of  specific  performance  parameters  for  arbitrary 
target  shapes  in  the  presence,  of  noise  and/or  clutter,  along  with  performance 
bounds  On  what  an  optimal  process  would  be.  Second,  there  is  a  means  to 
show  that  for  optimum  target  shape  extraction,  the  algorithms  must  account 
for  the  presence  of  clutter,  or  adapt  to  it. 
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An  overall  modeling  approach  for  image  processing  is  presented  in 
Section  2.  This  represents  the  approach  used  for  this  report  where  analytic 
results  were  obtained  for  a  simple  scene  and  model  as  in  reference  1, 
and  then  extended  by  increasing  the  complexity  of  the  model  while  maintain¬ 
ing  tractability  of  the  analysis. 

Section  3  provides  a  brief  background  of  the  techniques  used  from 
detection/estimation  theory  and  a  definition  of  concepts  used  in  defining 
performance  of,  e.g.,  an  object  segmentor. 

Sections  4-6  follow  the  modeling  philosophy  of  Section  2  by  extending 
the  basic  technique  of  computing  performance  bounds  to  more  complex 
models.  The  emphasis  is  on  improving  the  realism  of  the  model  while 
maintaining  the  tractability  of  the  analysis,  in  this  case,  the  tractability 
of  the  Fisher  information  matrix,  which  is  the  key  element  in  computing 
analytic  performance  bounds.  In  Section  4,  a  more  complex  target  model 
is  developed  and  expression  for  the  elements  of  the  Fisher  information 
matrix  derived  for  the  recommended  approach.  In  Section  5,  the 
complexity  of  the  scene  model  is  increased  by  introducing  clutter.  Two 
types  of  clutter  models  are  described  and  some  complex  drivations  are 
presented  to  show  that  analytic  results  can  be  obtained.  A  verification  of 
the  results  for  the  first  model  is  presented  by  showing  satisfactory  agree¬ 
ment  with  results  from  a  Monte  Carlo  computer  simulation.  In  Section  6, 
the  complexity  of  the  model  is  increased  by  introducing  a  sensor  model  that 
exhibits  detector  quantization  and  analytic  results  of  this  development  are 
presented. 

Section  7  contains  an  analysis  of  the  principal  features  of  the  ATAC 
segmentor.  Included  is  an  analysis  of  the  Sobel  edge  operator  with  an 
analytic  procedure  for  finding  the  probability  of  gaps  in  edges  and  an  experi¬ 
mental  verification;  and  an  analysis  of  the  selection  of  interior  points  to 
define  the  segmented  object. 

Recommendations  for  further  development  of  analytic  techniques  in 
image  processing  are  presented  in  Section  8. 

The  appendices  contain  some  of  the  more  detailed  and  complex  deri¬ 
vations  including  some  of  the  mathematical  and  statistical  manipulations 
and  concepts  that  have  been  found  useful  in  the  development  of  the  various 
analytic  results. 


2.0  IMAGE  PROCESSING  MODELING  PHILOSOPHY 


The  general  problem  of  modeling  an  image  processing  system  is 
illustrated  in  Figure  2-1  as  four  interrelated  tasks.  A  scene  is  sensed  by 
a  sensor  such  as  a  FLIR,  an  algorithmic  process  is  applied  to  the  sensed 
data  to  obtain  some  reduced  and  meaningful  data  and  an  objective  qualita¬ 
tive  measure  of  the  performance  is  obtained.  The  main  difficulties  in 
generating  the  performance  model  are  the  extreme  complexity  of  the  possi¬ 
ble  scenes  and  algorithms,  and  defining  appropriate  and  meaningful  per¬ 
formance  parameters.  In  order  to  obtain  a  tractable  performance  analysis, 
the  models  of  the  three  principal  components  (scene,  sensor,  process)  need 
be  simplified  to  account  for  the  particular  performance  parameters.  These 
depend  on  the  very  process  being  modeled,  which  in  turn  depends  on  the 
scenes  to  be  recognized. 

The  philosophy  for  modeling  the  ATAC  autocuer  used  in  this  effort 
was  to  start  with  a  simple  scene  model,  in  this  case  a  rectangular  target 
imbedded  in  noise.  Then  the  algorithmic  processor,  or  portions  thereof, 
are  examined  to  determine  its  basic  objective  and  which  kncv'n  analytic  tech¬ 
niques  (e.g.  detection/estimation  theory)  were  suitable  for  describing  its 
performance.  This  results  in  both  a  tractable  analysis  for  a  simple  case 
and  the  definition  of  appropriate  performance  parameters.  Then,  using  those 
performance  parameters,  the  complexity  of  the  scene,  sensor  and  process 
models  are  increased,  while  maintaining  the  tractability  of  the  analysis. 

This  modeling  approach  and  philosophy  does  not  result  in  attempting  to  model 
all  possible  real  world  scenes,  but  models  the  basic  scene  characteristics 
that  have  a  recognizable  effect  on  performance  of  the  algorithm  process. 
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MODEL 
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Figure  2-1.  General  Model  Diagram. 

This  process  is  iterated  into  more  complex  models  w  here  the  basic 
characteristics  are  modeled  w'ith  the  fewest  possible  parameters  and  in  a 
functional  form  that  maintains  a  tractable  performance  analysis. 


3.0  BASIC  ANALYTIC  TECHNIQUES 


The  basic  analytic  technique  used  in  the  analysis  of  the  simple 
rectangular  set  mentation  (1)  is  based  on  detection  and  estimation  theory  (2). 
This  approach  uses  the  concept  of  maximum  likelihood  estimation  to  derive 
performance  bounds  on  the  optinsum  segmentation  of  a  rectangle  imbedded 
in  white  noise  and  defines  the  best  possible  performance  against  which  any 
segmentation  algorithm  can  be  compared.  This  section  provides  a  basic- 
technical  summary  of  the  technique  for  both  sampled  and  continuous  wave 
fo  rms . 

In  order  to  develop  the  principles  of  maximum  likelihood,  it  is  con¬ 
venient  to  develop  the  principles  for  the  "classical"  estimation  problerr 
where  the  signals  presented  to  the  estimator  are  discrete  samples.  After 
developing  the  necessary  results  for  the  "classical"  estimation  problem 
these  principles  are  extended  to  the  case  where  the  observation  presented 
to  the  estimator  is  a  continuous  wacoform.  In  addition,  the  development 
will  be  for  the  case  where  the  observed  signals  are  one  dimensional,  since 
the  two  dimensional  case  follows  easily  from  this. 

The  use  of  the  concept  of  maximum  likelihood  is  significant  in  that 
it  provides  the  optimum  performance  of  the  detection/estimation  problem. 

With  a  knowledge  of  optimum  performance,  comparisons  with  specific 
algorithms  can  be  made  to  determine  if  significant  improvements  can  be 
made  and  possibly  show  where  or  in  which  direction  to  look  for  those  improve¬ 
ments.  In  addition,  there  are  powerful  results  in  detection/estimation 
theory  that  permit  computation  of  performance  bounds  for  these  optimum 
processes. 
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3.1  L'STIMATION  Tllt'ORY 


For  the  sampled  signal  case,  the  estimator  observes  the  N  discrete 
samples 

r(t.)  =  S(t.,  A)  +  n(t.)  (i  =  1.  2,  3,  .  .  .  ,  N)  (5-1) 

where  ^  is  tlie  feature  vector  representing  the  parameters  to  be  estimated; 
the  functional  form  of  S(t.,  _A)  is  known  but  the  K  x  1  vector  A  is 
unknown;  n(t^)  is  a  white  Gaussian  process  of  zero  mean  and  variance  N^/2. 
The  estimator  is  required  to  use  the  N  samples  to  estimate  the  pai-ameter  A. 

An  estimator  of  A  is  typically  denoted  by  A(r(t,  ),  rtt^),  .  .  .  ,  r(t-,)  I  A). 
It  is  incportant  to  note  that  the  estimator  is  a  function  of  the  observed  sam¬ 
ples,  although  we  might  write  the  estimator  simply  as  A. 

Noting  the  dependence  of  _A  on  the  N  observations,  which  are  random, 
it  is  clear  that  A  is  a  random  variable  with  an  expected  value.  If  this 
expectation  is  denoted 

^(A(r(t^),  ...,  .  .  .  ,  r(tj^)|  A)  =  A  (3-2) 

A  is  said  to  l)e  unbiased  wlien  this  equality  holds. 

Clearly,  it  is  reasonable  to  seek  an  estimator  which  will  give 
estimates  which  are  close  to  the  true  parameter  value.  Hence,  (whenever 
possiljle)  an  estimator  is  used  whose  expected  squared  error  around  any 
element  of  A  is  small;  that  is  (if  possible)  an  estimator  is  used  such  that 

/(a.  -  aH  is  small  for  each  i(i  =1,  2,  ...,  Kl  (3-3) 


whe  re 

q'  ^  .  T 

A  =  (a^,  ^2*  •••>  ~  ^2'  ^3*  •••>  * 

If  an  unbiased  estimator  is  used,  the  well  known  Cramer-Rao 
(Ref.  2)  inequality  can  be  used  to  give  an  iclea  of  how  much  squared  error 
there  will  be  in  using  the  "best  "  unbiased  estimatcir. 


In  order  to  apply  the  Cramer-Rao  inequality,  first  define  the  Fisher 
information  matrix,  J,  which  has  elements 


J.. 

ij 


/  din  P(r(t^),  ...»  r(tj^)) 

\ 


din  P  |r{t2 ,  . . .  , 


da. 

J 


(3-4) 


where  P(r(tj^),  ^■(t^),  ...,  is  the  probability  density  or  likelihood  of 

the  N  observations.  The  Cramer-Rao  inequality  is  then  seen  to  be  (2) 


f  ((a.  -  =  Var  (a.  -  a.)  >  (3-5) 

where  is  the  ii*"^  element  of  the  KXK  matrix  J  An  unbiased  estimator, 

A 

A,  is  said  to  bo  efficient  if  equality  holds  in  (3-5).  Because  of  the  value  of 
the  Cramer-Rao  inequality  in  determining  performance  bounds,  the  Fisher 
information  matrix  becomes  an  important  concept  in  the  analysis  of  the 
remaining  sections. 

It  must  be  pointed  out  that  there  is  no  reason  to  assume  that  an  effi¬ 
cient  estimator  will  always  exist,  and  in  fact,  it  is  easy  to  find  examples  for 
which  an  efficient  estimator  does  not  exist.  However,  if  an  efficient  esti¬ 
mator  A  exists,  then  (with  probability  of  one)  it  will  be  given  by  the  solution 
to 

f  (A)  =  max  P(r(t  J  ),  ...,  r(tj^)|  ,  (3-6) 

which  is  the  maximum  likelihood  estimate. 

Since  maximum  likelihood  estimators  maximizes  the  likelihood  of  the 
observed  samples,  it  is  desirable  to  use  them  whether  an  efficient  estimator 
exists  or  not.  Then,  the  estimates  will  in  general  be  quite  good  since  for 
certain  general  conditions  it  has  been  shown  that  maximum  likelihood 
estimates 

(1)  converge  to  the  correct  parameter  values  as  N  ->  33 

(2)  become  efficient  as  N  -*■  “> 
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Using(3-1),  (3-6)  and  the  assumptions  about  M(t^)  the  likelihood  is  seen  to  be 


i{A)  =  P(r(tj),  ....  r{tj^))  =  n  P.(r(t.)|A) 

j=l 


(3-7) 


where  the  Pj{r(tj)|  A)'s  are  the  probability  densities  of  the  r(tj)'s.  Then 

where  is  the  noise  spectral  power.  It  is  usuadly  easier  to  find  the  value 
of  A  which  maximizes  In  in  £(A)  rather  than  t{A).  Therefore,  it  is  usually 
best  to  work  with  the  log-likelihood. 


Aj(A)  =  ln(A) 


(3-9) 


which  reduces  to 


N  1  2  2  '^  / 


(3-10) 


Note  that  in  (3-10)  the  first  two  terms  do  not  depend  on  the  parameters  to 
be  estimated.  Consequently,  the  first  two  terms  would  typically  be  dropped, 
and  the  log -likelihood  would  be  redefined  to  be 


^2^-^  "  No 


N  / 

S(t.,  A)  |r(t.)  - 
i=l  \ 


S(t.,  A) 


(3-11) 
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Equation  (3  -11)  is  the  form  of  the  log -likelihood  that  would  typically 
be  used  in  estim.ating  A  and  in  computing  the  Fisher  information  matrix. 
For  the  continuous  wave  form  case,  this  equation  converges  to 


A 


,  _ 

"  No 


r(t)  S(t, 


A)  dt 


No 


S^(t,  A)  dt 


(3-12) 


This  is  the  form  of  the  log -likelihood  which  would  be  used  to  compute 
the  maximum  likelihood  estimator,  the  Fisher  information  matrix,  or  any 
other  quantity  associated  with  maximum  likelihood  estimation  in  the  con¬ 
tinuous  case. 

3.2  PROBABILITY  OF  ANOMALY 

The  do\-clcpmenhs  in  Section  3.  1  provide  the  basic  means  for  develop¬ 
ing  tight  performance  bounds  when  the  process  results  in  output  near  the 
true  results.  For  example,  in  the  problem  of  segmenting  an  object  from 
the  background  those  developments  can  be  used  to  obtain  good  performance 
bounds  w'hen  that  object  is  segmented,  but  not  exactly.  To  describe  com¬ 
plete  performance  it  is  necessary  to  define  an  additional  criteiion,  termed 
the  probability  of  anomaly.  This  is  the  probability  of  detecting  falsely. 

In  the  segmentation  example,  it  is  the  probability  that  the  true  object  is 
not  formed,  but  rather  the  scgmcntor  "finds"  something  in  another  part  of 
the  seem;  due,  e.g,,  to  noise.  This  situation  occu rs  when  the  likelihood 
function  in  an  area  containing  only  noise  exceeds  the  likeliliood  turn  lion 
of  the  area  containing  the  object  itself. 

The  com.putation  of  this  probability  of  anomaly  uses  the  same  basic 
concepts  as  the  performance  bounds,  that  is,  the  likelihood  equation  |3ro\-ides 
the  ki'y  parameters  used  in  the  analytic  computations.  Examples  of  this 
process  are  gi\'en  in  later  sections. 
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4.0  TARGET  MODELS 


The  first  problern  in  the  analysis  of  target  classification  algorithrns 
using  the  techniques  of  maximum  likelihood  estimation  is  to  obtain  simple 
target  models  that  are  tractable  analytically.  By  this  it  is  meant  that  the 
Fisher  information  matrix  (Section  3) 

T  -  A.  ff  y>  i) 

'ij  -  No  JJ  8f. 

can  he  easily  inverted  wliere  I(x,  y,  J)  is  the  intensity  function  of  the  target 
model,  and  that  the  feature  vector  _f  describing  the  target  has  the  smallest 
dimension  possible. 

4.1  APPROACHES  EXAMINED 

In  this  section  a  number  of  possible  representations  of  I(x.  y,  f)  are 
described. 

Z-D  Legendre  Polynomial  Representation 

These  polynomials,  P(  ■  ),  are  orthogonal  on  a  square  window  and  in 
principal  can  be  used  to  approximate  any  continuous  shape.  If  the  target 
intensities  are  given  by  l(x,  y)  then  the  least  square  approximation  by  the 
polynomials  is  given  by  minimizing  the  expression 

2 

i  j 
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for  each  i  and  j.  Since  these  polynomials  are  orthogonal,  the  minimum  is 
given  by  the  coefficients  a. .  defined  by 


l(x,  y)  P^(x)  Pj(y)  dxdy 


(4-3) 


Since  the  polynomials  P^(x)  and  Pj(x)  are  polynomials  in  x  and  y  of  degree  i 
and  j,  these  coefficients  are  just  linear  combinations  of  the  various  moments 
of  the  target  intensities.  The  motivation  for  this  approach  was  that  many 
target  classification  approaches,  including  the  ATAC  FLIR  autocuer  use 
moments  as  target  features  in  the  classification  process.  In  order  to  deter¬ 
mine  how  well  these  polynomial  moments  can  approximate  the  shape  of  a 
target,  an  experiment  was  performed  on  a  typical  tank  shape.  The  polynomial 
moments  a.,  were  computed  using  Equation  (4-3)  and  the  corresponding 
shape  approximation  reconstructed  using  the  summation  term  in  Equation 
(4-2),  for  several  cases  of  different  orders  of  moments.  The  results  arc 
shown  in  Figure  4-1,  where  it  is  apparent  that  a  good  approximation  can 
require  moments  to  the  20th  order.  This  would  require  over  200  coefficients 
for  target  representation,  which  is  too  high  to  be  practical.  In  addition,  the 
coefficients  are  not  immediately  scale  and  rotation  invariant. 

The  principal  advantage  of  this  approach  is  that  due  to  the  ortho¬ 
gonality  of  the  polynomials,  tlie  Fisher  information  matrix  (4-1)  becomes 
diagonal  and  need  not  be  inverted.  However,  this  did  not  prove  sufficiently 
advantageous  to  adopt  tliis  approach. 

Composite  Rectangle  Model 

The  second  model  considered  was  the  composite  rectangle  model 
(Figure  4-2).  In  this  approach  an  arbitrary  target  shape  is  approximated 
by  rectangales  of  various  siices.  The  advantage  of  this  approach  is  that 
any  shape  can  be  reasonably  approximated.  The  disadvantage  is  that  it  is 
not  tractable  mathematically.  This  is  because  each  rectangle  must  be 
adjacent  to  at  least  one  other,  causing  many  non-diagonal  elements  in  the 
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Figure  4-1 .  Reconstruction  using  legendre  moments. 


Figure  4-2.  Composite  rectangles. 


Fisher  matrix.  It  also  is  difficult  to  model  the  adjacency  of  rectangle 
edges  mathematically  for  cirbitrary  configurations. 

Composite  Equal-Angled  Triangular  Model 

This  type  of  model,  as  illustrated  in  Figure  4-3,  can  approximate 
many  shapes  and  is  mathematically  more  tractable  than  the  composite 
rectangular  model  since  triangles  have  a  more  easily  described  relationship 
to  each  other  than  do  rectangles.  The  disadvantage  is  that  sonie  shapes, 
those  for  which  the  radius  vector  is  not  a  single -valued  function  of  angle, 
cannot  be  represented  with  this  model. 
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Figure  4-3.  Composite  equal-angled 
triangular  model. 


General  PolyKOt^al  Model 

This  approach  is  to  approximate  any  shape  by  specifying  points  on 
the  boundary  and  attaching  them  with  straight  line  segments.  An  example 
of  this  approach  is  shown  in  Figure  4-4.  Any  shape  can  be  ap)^ roxiniated 
in  this  way  realistically  using  relatively  few  parameters.  In  addition,  this 
approach  exactly  describes  a  rectangle  so  it  is  a  direct  extension  of  rec¬ 
tangular  targets  to  more  complex  shapes.  These  features,  as  well  as  the 
Lractability  ot  this  representation  resulted  in  the  selection  of  this  representa 
lion . 

The  results  for  the  general  polygonal  case  are  as  follows.  Let  the 
Fisher  matrix  (4-1)  be  applied  to  the  case  of  a  polygon  of  N  sides  whose 
vertices  are  specified  by  coordinates  (x^, 
shown  in  Figure  4-4. 

It  is  assunied  that  inside  the  polygon  I{x,  y,  _1)  =  AT,  while  outside 
Uu;  polygon  1  =  0,  except  at  points  nearer  to  the  boundary  than  a  small  dis¬ 
tance  ±  h/Z  (assumed  to  be  small  compared  with  any  side  of  the  polygon). 

It  is  assunied  that  within  this  small  region  of  width  1  i-aries  linearly  from 
its  value  AT  inside  the  boundary  to  its  value  zero  outside  the  boundary. 


{^2-  '^2* 
(X,.Y,) 

(Xq.  Vq)  =  (X^.  Y^) 


Figure  4-4.  General  polygon  model. 


It  is  further  assumed  that  the  first  shape  feature  =  AT,  and  that 
for  all  parameters  f^,  L,  i,  j  >  1,  I  depends  upon  ^  only  through  the  vertices 
(x^,  y^)  =  1.  2,  ....  N.  That  is 


^n  ^n^^2’  ^3’  *  ‘  ‘ 


''^n  ~  '^n^^Z*  ^3’  •  •  •  ’ 


M  <  2N  +  1 


Aside  from  these  restrictions,  the  f^'s  can  be  chosen  arbitrarily  so 
long  as  they  are  independent. 

For  example,  two  coordinates  might  specify  a  centroid,  or  other 
fixed  point,  to  which  all  other  rotations  or  forms  are  to  be  referred,  one 
coordinate  might  specify  rotation  in  the  plane,  two  coordinates  could  refer 
to  aspect  angles  of  the  three  dimensional  object  being  projected  onto  the 
image  plane . 

Then,  to  a  good  approximation  (see  Appendix  G) 
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where  are  coordinates  of  the  midpoint  of  the  nth  side 


X  +  X  , 

—  n  n-1 

X  =  - :z - 

n  2 


y.  = 


y  +  y  1 

’ n  ^n-l 


n  2 


(Ax^,  are  the  relativ'e  coordinates  of  the  nth  side,  1-^  is  the  length 

of  the  nth  side, 

1  /-> 


L  = 
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(Ax  +  (Ay 
n  ^n 


0^  is  the  angle  (with  respect  to  the  positive  X  axis  made  by  the 


nth  side, 


e  =  'an 
n 


1  ^y 

-1  ^n 


Ax 


sin  0  = 


n  L 


Ax 

cos  0  =  — — 

n  j.. 


for  i,  j  >1,  the  elements  of  the  Fisher  information  matrix  are 
,  N 
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The  quantities  can  also  be  written  in  the  forn. 
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n  8f . 


=  L 


dAy  dAx 

n  .  ,  n 

cos  U  — — —  -  sin  0 - 
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It  has  been  verified  that  these  results  reduce  t  orrectly  to  those 
given  in  reference  1  for  the  case  of  a  rectangle. 

4.2  SIGNIFICANCE  OF  THESE  RESULTS 

These  results  extend  tlie  past  work  to  the  case  of  an  arbitrary  polygon 
depending  upon  an  arbitrary  set  of  feature  parameters.  The  only  requirement 
is  that  the  vertices  of  the  polygon  be  known  differentiable  functions  of  the  fea¬ 
ture  parameters.  Tlie  two-dimensional  integral  defining  the  Fisher  informa¬ 
tion  matrix  (FIM)  has  been  converted  into  a  relatively  simple  one -dimensional 
sum  over  the  sides  and  vertices  of  the  polygon. 

Algebraic  oi-  numerical  inversion  of  the  FIM  can  be  obtained  for  any 
case  of  practical  interest  to  obtain  bounds  on  the  variances  of  the  feature 
parameters. 

Tile  fact  that  the  model  is  an  arbitrary  N -sided  polygon  makes  it 
applicable  to  a  wide  variety  of  military  targets  and  clutter  objects.  As  loJig 
as  targets  and  clutter  objects  don't  overlap,  the  FIM  will  be  I'cadily  scii^ra- 
ble  into  corresponding  constituent  parts.  For  example,  an  actual  silouette 
of  a  target  of  interest  could  be  represented  by  selecting  points  on  tlie 
boundary  that  give  an  arbitrarily  accurate  representation  of  the  sliape.  , 

The  performance  bounds  for  optimum  segmentation  for  this  representation 
can  then  be  computed  using  the  results  of  this  section  applied  to  the  coordi¬ 
nates  of  these  points.  These  results  can  be  compared  for  differing  shape 
approximations  to  find  the  simplest  representation  of  the  target  that  does 
not  significantly  alter  the  analytic  results,  yielding  an  acceptable  target 
model. 
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5.0  CLUTTER  MODELING 


This  section  presents  analysis  and  results  on  modeling  the  effects 
of  simple  clutter  on  target  segmentation  p  roces ses .  Sections  5,1  and  5.2 
present  two  clutter  models,  Section  5.3  shows  results  of  the  presence  of 
clutter  on  target  segmentors.  Section  5,4  gives  results  of  a  Monte  Carlo 
simulation  to  %’erify  the  analytic  bounds  on  piobability  of  anomaly  (derived 
in  Appendix  11)  and  Section  5.5  shows  an  example  of  what  probabilistic 
textured  clutter  could  look  like. 

As  a  basic  assumption  it  is  assumed  that  clutter  is  any  object  in 
tile  I'OV  thru  is  not  a  target  or  noise.  The  only  restriction  placed  on  the 
clutter  is  tliat  it  never  masks  the  target  although  target  masking  of  clutter 
is  allovK  ed,  and  that  if  there  are  multiple  clutter  objects  in  the  FOV  tliey 
nes'cr  overlap.  For  any  arbitrary  shaped  clutter  object  in  the  plane  the 
clutter  object's  shape  can  be  approximated  by  the  polygon  method  described 
in  the  previous  section.  That  is,  one  could  connect  vertices  by  straight 
lines  to  approximate  the  shape  of  the  clutter  object.  The  vertex  locations 
can  be  assumed  to  be  unknown  or  random.  More  detailed  assumptions  will 
be  made  later,  in  order  to  get  a  realistic  but  mathematically  tractable 
clutter  model. 

5.1  CLUTTER  MODEL  1 

In  addition  to  tlie  tai  gct  parameters  that  must  be  estimated,  assume 
that  tile  vertices  of  each  clutter  oijject  arc  unknown  and  must  be  estimated 
and  that  the  intensity  offset  of  each  clutter  object  is  unknown  and  must  also 
be  estimated.  The  numijer  K  of  clutter  objects  is  assumed  known. 
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Figure  5-1  illustrates  the  setup  for  K  =  3.  This  can  oe  represented 


(x,  y)  =  AT  S{x,  y)*  P(x,  y) 

/  K 


^  Pj  (l  -  S(x,  y))  C^(x,  y)  >:=  P(x,  y)  +  n(x,  y) 


=  observed  intensity  dis'ribution 
,  (  1  if  X,  y  is  on  target 


S(x,  y)  = 


I  0  otherwise 

J  1  if  X,  y  is  on  clutter  object  j 


0  otherwise 


Figure  5-1.  One  target,  three  clutter  objects  and  white  noise. 


AT  =  target  intensity  offset 


nix,  y) 


intensity  offset  of  j  clutter  oliject, 
zero  mean  Gaussiaui  white  noise  process, 


Then  in  this  case  S^(k,y)  =  S(x,y)  and  C^(x,y)  »  Cj(x,y)  and  (5-3)  reduces  to 

A'  =  JJ  AT  S(x,  y)  |l(x,  y)  -  dxdy 

°  FOV 

+  ^  ^  JJ |(JLj  (l  -S(x,  y))  C.(x,  y)j  |l(x,  y)  dxdy 


(5-4) 

Note  that  this  is  tlie  sum  o£  K+1  matched  filters,  one  matched  to  the  target 
and  K  matched  to  tlie  K  clutter  objects. 

Due  to  the  separability  of  the  log  likelihood  function  into  target  and 
clutter  components  it  can  be  shown  that 


J  =  (J..)  = 
D 


(5-5) 


where  is  the  PxP  square  matrix  corresponding  to  a  p  parameter  target 
and  C  is  a  QK  xQK  square  matrix  corresponding  to  KQ  parameter  clutter 
object. 

For  a  rectangle  target  as  in  (1)  and  the  target  parameters  to  be 
estimated  are  as  in  (1),  it  is  not  hard  to  show  that  is  the  5x5  Fisher 
information  matrix  computed  in  reference  1,  and  C  is  the  inform.ation 
matrix  of  the  clutter.  Since 


(5-6) 


n  C 


the  same  target  estimation  bounds  apply  here  as  in  (1). 
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5.2  CLUTTER  MODEL  11 

This  clutter  model  is  similar  to  model  I.  All  the  assumptions 
of  model  I  apply;  however,  model  I  is  extended  so  that  each  clutter  object’s 
offset  intensity  is  random.  That  is. 


K 

I(x,  y)  =  AT  S(x,  y)  +  (1  -  S(x,y))  Cj(x,  y)  p 

j=l 


K 

r  (1  -  S{x,  y))  Cj(x,  y)  T^  +  n(x,  y) 

j=l 


(5-7) 


where  all  variables  are  the  same  as  in  model  I  except  that  Tj's,  which  are 
assumed  to  be  independent  Gaussian  with 


^(T.)  =  0,  Var{T.)  =  <7^,3=  1,  2,  3, 
For  this  model  the  total  noise  is  non-white 


K. 


K 

y)  =  (1  -  S(x,  y))  Cj{x,y)  Tj  +  n(x,  y) 

j=l 


The  autocorrelation  of  the  non-white  noise  is 


S'{ri^(x,y)  n^(u,v)|  =  K^(x,  y,u,v) 

K 

=  <7^  (1  -  S{x,y))  Cj(x,y)  (5-8) 

i=l 

•  (1  -  S(u,v))  Cj{u,v)  +  —  6(x-u,  y-v) 
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Let 


Kg(x,  y,  u,  v) 


K 

0-^  (1  -  S{x,y))  C^(x,y)  •  (1  -  S(u,  v))  C^{n.  v) 

j=l 

(5-9; 


K  (x,  y,  u,  v)  is  a  separable  kernel  as  in  reference  3  (Section  7,3). 
Since 
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FOV 


(1  -  S(x,  y))  ^ 


C.(x.v)  C,(x,y) 


dxdy  = 
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tH 

where  A.  is  the  area  of  the  i  clutter  object,  set 
1 

<J>^(x,  y)  =  (1  -  S(x,  y))  Cj(x,y)/^AT  ,  j=l,  2,  3,  .  .  , ,  K 

X.  =  0-^  A.  ,  j=l,  2,  3,  . .  .,K 

) 

X.  -  0  ,  j>K 


and 


^  A  is  area  of  target. 

^K+L  ’  s 


(5-lOa) 


(5 -10b) 
(5-lOc) 
(5-lOd) 


(5-lOe) 


Let 
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m(x,  y)  =  -(at  S(x,  v)  + 

j=l 


p(l  -  S(x,  y)) 


Cj(x,y); 


(5-11) 
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Then  the  log -likelihood  is  (see  reference  3,  Chapt.  6) 
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i=l  ^  i=l 


No/2 


i=K+l  -i=l  i=l  1 


.  +  No/2 


i  V 

'  2  Z_^  No/2 

L=K+1 


(5-12) 


where 


r.  = 
1 


II 

FOV 


C.(x,y) 

I(x,  y)  (1-S{x.  y))  - dxdy  i=  1 ,  2,  3,  .  .  . ,  K 

yjA. 


and 


K+1 


/ /  I(x,  y)  dxdy  , 

4r 


FOV 


ff  C(x,y) 

j=  JJ  rn(x,  y)  (1  -  S(x,  y))  - ; -  dxdy 


FOV 


=  \/a. 


i  —  \  K 
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Then  the  log -likelihood  (5-12)  reduces  to 


K  2  .  K 

=  -L  V  i  2  M 

‘  .^a..no/z  J 


No/2 


r . 
J 


FOV 


iZ-  (-^)- 


K 


i=l 


AT^  A 


iE 


j  _  a  A  j  +  N  o  /  2 


No 


To  conipute  the  Fisher  information  matrix,  set 


K  ‘  cr'-  A. 


h(x,y,u,v)=  C^Cx.y) 


i=l  \  2 


(1  -  S(u,  v))  C^(u,  v)  , 


Q^(x,y,u,vl  =  (6(x  -  u,  y  -  v)  -  h{x,  y,u,v)) 


and 


K 


K^(x,y,u,v)  0-^  A.  (1-S(x,y))  C.(x,y)  (1-S(u,v))  C.(u,v 


i=l 
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(5-12b) 


(5-13) 

(5-14) 

(5-15) 


Then  from  Reference  3,  Chapter  7,  the  elements  of  the  Fisher  information 
matrix  are 


kf 


//// 

FOV  FOV  L 


22^  .  C!>,y,u,v)  • 


9f  { 


I  ak^{x,  y,  u,  v)  aQ^(x,  y,  u,  v) 


af, 


8  f  f 


dxdy  dudv' 


(5-16) 


In  this  section  it  was  shown  that  for  this  model  there  is  a  closed  form 
expression  for  the  likeliliood  (equation  5-12b),  which  is  needed  to  compute 
probability  of  anomaly.  Also,  there  is  an  expression  for  the  elements  of  the 
Fisher  information  matrix  (5-16).  Applications  of  these  formulae  to  compute 
the  probability  of  anomaly  for  scenes  with  simple  clutter  can  be  found  in 
Appendix  11. 

5.  3  EFFECTS  OJ  CLUTTER  ON  TARGET  SEGMENTATION 

In  order  to  evaluate  the  effects  of  clutter  on  target  segmentation,  it 
will  be  assumed  that  there  is  a  target  and  a  clutter  object  of  known  shape 
but  unknown  size  and  location.  With  this  assumption  it  is  possible  to 
determine  probability  of  anomaly  for  three  cases  of  interest.  (See 
Ajjpendix  11.) 

1.  There  is  a  target  and  a  clutter  object  in  the  scene,  and  the 
segmentor  is  optimized  for  this  situation. 

2.  There  is  a  target  and  a  clutter  object  in  the  scene,  but  the 
sementor  is  optina^ed  for  a  target  only  in  the  scene. 

3.  There  is  only  a  target  in  the  scene  but  the  segmentor  is  opti¬ 
mized  for  target  and  clutter  object  in  the  scene. 

Figure  5-2  shows  the  probability  of  anomaly  for  these  cases,  where 
the  target  is  a  5  x  5  square  and  the  clutter  is  a  5  x  10  rectangle,  as  a  func¬ 
tion  of  signal  to  noise  ratio.  These  curves  show  clearly  that  the  segmentor 
should  be  constructed  to  account  for  the  presence  of  clutter  to  get  the  best 
performance.  That  is,  any  segmentation  algorithm  shovild  specifically 
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CNR  -  1 

target  is  6  X  6  SQUARE 
CLUTTER  IS  6  X  10  RECTANGLE 

TARGET  AND  CLUTTER 
IN  FOV  SEGMENTATOR 
OPTIMIZED  FOR  TARGET 
AND  CLUTTER  IN  FOV 

TARGET  AND  CLUTTER 
IN  FOV  SEGMENTATION 
OPTIMIZED  FOR  TARGET 
IN  FOV 

TARGET  IN  FOV  SEGMENTA¬ 
TOR  OPTIMIZED  FOR  target 
AND  CLUTTER  IN  FOV 


Figure  5-2.  Effects  of  clutter  on  probability  of  anomaly. 

account  for  or  adopt  to  the  presence  of  clutter  objects  in  the  vicinity  of 
targets  in  order  to  achieve  good  performance.  Otherwise,  if  there  is  a 
clutter  object,  performance  will  be  limited  to  be  no  better  than  the  second 
curve  in  Figure  5-2, 

5.4  VERIFICATION  OF  PROBABILITY  OF  ANOMALY  BOUNDS 

As  stated  above,  the  derivation  of  analytic  bounds  on  probability 
of  anomaly  for  segmentation  of  a  target  in  the  presence  of  clutter  is  given 
in  Appendix  H.  Some  verification  of  the  bounds  using  Monte  Carlo  estima¬ 
tion  of  probability  of  anomaly  is  desirable  to  verify  the  analytic  bounds. 
Since  a  general  verification  of  the  principal  result  (Equation  H26)  was  not 
feasible  due  to  the  Ingli  dimensionality  of  the  problem  and  the  consequent 
high  computer  processing  cost,  a  modified  version  of  the  general  problem 
was  selected  for  verification.  The  basic  simplifying  assumptions 
were; 

1.  The  target  is  known  to  be  a  square. 

2.  The  clutter  object  is  known  to  be  a  rectangle  of  width  and 
height  2 

3.  The  vertical  (in  the  scene)  location  of  the  target  and  clutter 
is  known. 

These  assumptions  reduced  the  dimensionality  of  the  process,  allowing 
for  a  reasonable  verification  to  be  performed.  The  results  of  this 


© 

© 
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verification  are  shown  in  Table  5-1  where  it  is  seen  that  results  from  the 
analytic  bound  computation  and  Monte  Carlo  Simulation  are  similar. 


TABLE  5-1.  COMPARISON  OF  ANALYTIC  BOUND  ON  PROBABILI  TY 
OF  ANOMALY  WITH  MONTE  CARLO  SIMULATION 


Signal  to  Noise 
Ratio 

Analytic 

Bound* 

”  - ■  “  1 

No.  of  Anomalies 
from  Simulation 

No.  of 

T  rials 

.  1 

1 

10 

10 

.2 

0.402 

i, 

10 

.  3 

0.  35 

4 

10 

,  5 

0.254 

2 

10 

1 

0.  013 

0 

30 

1 

7  x  10""^ 

0 

10 

4 

0 

0 

10 

»:-CNR  =  1  in  Eq.  (H26) 

5.5  TEXTURED  CLUTTER 

The  results  of  the  previous  two  sections  were  for  the  first  type  of 
clutter  objects  in  the  target  scene;  a  natural  extension  of  this  analysis  would 
include  consideration  of  tlie  textured  background  tyue  of  cluttei-.  An  initial 
approach  would  be  to  assume  that  the  textured  background  can  be  repre¬ 
sented  by  a  wide-sense  stationary  random  field.  In  this  case,  the  textured 
background  can  be  represented  by  the  equation. 


nj3  (x,  y)  =  g(x  -  1,  y)  +  o  g  (x,  y  -  1) 
a  ‘  ^ 


-  Uy  g(x  -  1,  y  -  1) 


(5-17) 
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where  a 

X 


by 


and  a 

V 


are  the  texture  features  parameters 


and  g{x,  y)  is  given 


g(x,y)  =  S(x,y)  +  C(x,  y) 


(5-18) 


+  (1  -  C(x,  y))  (1  ~  S(x,y))  (x,y)  +  n(x,y) 

a 

The  remaining  terms  of  the  latter  equation  are  as  defined  in  Appendix  H, 

An  example  of  this  t^'pe  of  textured  background  is  shown  in  Figure  5-3, 
as  generated  by  computer  simulation  of  the  above  equations  with  no  target 
or  clutter  objects. 

An  extension  of  this  concept  of  clutter  to  verify  its  realism  to  real 
Fl.lR  image  effects  and  a  possible  extension  of  th.e  analytic  results  to  this 
case  is  a  futiire  research  task. 


es  SO  rs  too 

JUlilll  111  ■  H  I  I  ■  ■  I  1 1 


Figure  5-3.  Example  of  Markov  random  field. 


6.0  SENSOR  MODEl, 


The  analytic  results  obtained  in  the  previous  sections  did  not  include 
a  realistic  sensor  model,  in  particular,  the  effects  of  sensor  spatial  quanti¬ 
zation.  In  this  sectioji,  the  basic  analytic  teclmiques  for  computing  per¬ 
formance  bounds  for  an  optimum  segmentation  of  a  simple  scene  are  extended 
to  account  for  this  improved  sensor  model. 

A  natural  approach  is  to  compute  all  pertinent  quantities  (Eishers 
information  matrix,  likelihood  equation,  performance  measures,  etc.)  in 
terms  of  the  spatially  sampled  intensity'  distribution  which  is  present  to  the 
object  segmenter.  Wo  could  then  vary  the  sampling  rate  and  obsei  ve  tlie 
resulting  behavior  of  any  performance  measures  being  studied. 

In  order  to  apply  this  method  to  the  model  presented  in  Reference  1, 
first  the  statistical  clia  racte  ristic  s  of  the  sampled  intensity  distribution 
IS  determined;  then  the  likelihood  equation  is  derived,  from  which 
CranuT-Uao  anti  probability  of  anomaly  Ijountls  can  bi-  determined. 

The  tlerivation  of  tile  statistical  characteristics  of  the  sampleti 
intensity  distribution  is  as  follows: 

Let  l(x,y)  be  the  analog  int<-nsity  distribution  on  the  image  plane, 

that  IS  , 


Ux,  y  )  s  S(n,  y)  -t  n(x,  y) 


(h-1) 


win: re 


f 


S(x,  >  )  =  AT  •  Sj(x)  •  S  ,(y) 


(<■-2) 


1 


is  the  target,  and 


Sj(x)  - 


(o 


I  X  -  c  j  <  L  /2 
I  X I  x' 

otherwise 


j  1  |y  -  Cyj  1  Ly/2 
|o  otherwise 


(6-3) 


(6-4) 


and  n(x,  y)  is  a  zero  mean  Gaussian  white  noise  process  with  spectrum 
amplitude  . 

Assuming  a  detector  of  dimension  r  x  r  and  of  uniform  responsivity, 
the  signal  available  to  the  object  segmentor  from  a  detector  whose  center 
is  at  (x,y)  is  proportional  to 


whe  re 


Ij(x,  y) 


l(x  -  u,  y  -  v)  P(u,  v)  dudv 


P(u,v) 


|u|  <  r/2,  |v|  <  r/2 
othe  rwise 


(6-5) 


(6-6) 


r  r  T  . 

Assuming  that  tlie  samples  are  taken  at  AX  =  A  +  -^  a-nd  AY  =  A  +  ~  units 
apart,  the  sampled  image  that  the  object  segmentor  will  process  is;  for 
1  <  1,  j  i  N:  I^(x.,yj)  Ij  iAx  (i  -  1  2),  Ay  (j  -  1  2)) 


r/2  v/Z 

S(Ax(i  -  1/2)  -  u,  Ay(j  -  1/2))  •  P(u,  v  )  dudv 

•^-r/2 


r/2  Vj 

L  f 


j/^  ];/2 


I  I  P(u,v)  n(  Ax(i  -  1/2)  -  u,  Ay(j  -  1/2)  -  v)  dudv 

•'r/2  •'r/2 


(6-8) 


t 


t.-2 


The  mecLn  of  the  sampled  image  is 


u-j  =  ^  (is(^.yj)) 


(6-9) 


whicli  can  be  reduced  to 


AT 

u^.  =  — ^  Area  (Pixel  (i,j)  in//') 

where  is  the  set  of  all  points  that  lie  on  the  target.  The  variance  of 
the  sampled  image  is 


(6-10) 


r/2  r/2 

Var  I^(x.,y^)  =  Var  /  /  (P(u,  v)  n(Ax(  i  -  1/2)-  u  , 

•'-r/2  •'-r/2 


Ay(j  -  1/2)  -  v))dudv 

I 


■■  /  / 

2  -'-r/2 


P(\i,  v)  n{Ax{i  -  1/2)  -  u. 


Ay(j  -  l/2)  -  v)  dudv 


^valuation  of  the  integrals  results  ii 


(6-11) 


N 


Var  I^(x.,y_.)  = 


2r 


(6-12) 


Having  coni'pletely  specified  t'ne  statistii.  al  <.  ha  ra cte  ri stic  s  of  the 

sampled  image  1  (x.,y.),  this  segmentation  problem  can  be  treated  almost 
s  1  j 

exactly  as  the  continuous  segmentation  problem  in  Reference  1.  The  object 

segmenter  observes  Ig(x-»  Yj)  i,  j  =  1,  2,  3,  ...  ,N  and  is  required  to  estimate 

C  ,  C  ,  L  ,  L  and  AT. 
x'  y’  x’  y 
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As  in  Reference  1,  the  Fisher  information  matrix  is  used  to  obtain 
estimation  error  bounds.  First,  however,  the  likelihood  equation  must  be 
determined  so  that  Fisher's  information  matrix  can  be  obtained. 

The  likelihood  is; 


e  = 


(I  (x.,  y. )  -  u 
s'  i’ 


Taking  the  log,  the  log -likelihood  is  seen  to  be 


(6-1  i) 


r'^N, 


ij 


-  2  l^(x^,y^)  |i 


(6-14) 


Noting  that  the  first  two  terms  are  independent  of  the  parameters  to  be 
estimated,  we  write  the  1  og-liklihood  as 


Ai  = 


r  N, 


N,,r 


-Z 


'  N 

N 

E 

E 

Z 

^  1  =  1 

j  =  i 

'  N 

N 

E 

E 

1  =  1 

.i  =  i 

N 

N 

|x. .  -  Z  I  (x.,  y. )  u. 

s'  i'  •’j' 


(6-15) 


(6-16) 


Area  (Pixel  (i,  j)on'y')  Ig(x.,yj) 


.^1  ,  =  1  " 


(6-17) 


/AT 

Ar^ 


Area  (1  'ixel  (i,  j)  on 
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Using  (6-15)  -  (6-17)  to  compute  the  Fisher  information  matrix, 
which  by  definition  is  given  by 


J  =  (Jij)  =  <' 


IdAi  aAi\ 

\af.  '  9f.  / 

j  ' 


(6-18) 


where  Aj^  is  the  log -likelihood  and  f^,  F  are  parameters  to  Ije  estimated, 
results  in  the  expression 


N  N 


•  1  _  1  ' 


8|jl.  . 

) 

’  9f,. 


9|ii- 


i=l  j=l 

and  after  some  manipulation, 


(6-19) 


I  9A 


i=l  j=l 


(6-20) 


or, 


N  N  /at  »  Area  (i Uxel  (i,  j)  on // j\ 

■"  •" '  =  (if)  E  E  ^ f - - 


%  ’ 


i=l  j=l 


I  AT  •  Ar 


Area  (Pixel  ii,  j)  on 


"■) 


9f 


e 
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in  order  to  obtain  closed  form  results  the  following  assumptions  will 
be  made: 


(i) 


Area  (Pixel  {i,j)on//')  1 

2  “2 


(6-22) 
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for  any  pixel  which  is  on  the  target  boundary, 

(ii)  The  number  of  pixels  on  the  target 


L  L 

N  4  ---y  ^  .  (6-23) 

(r  +  A)^ 


(iii)  The  number  of  pixels  in  the  target  interior 


(iv)  The  number  of  pixels  on  the  target  boundary 


(6-24) 


or 


N^^Nt  -Ni 


/  2  \ 

Ir  ,  \ 

(r  +  a) 

K  +  Ly  -  2  (r  +  A)) 

(6-25) 


(6-26) 


(v)  The  number  of  pixels  on  the  target  boundary  parallel  to  the 
X-axis 


Ng  4  2Ly(r  +  A) 

X 


(6-27) 


(\-i)  The  number  of  pixels  on  the  target  boundary  parallel  to  the 
Y -axis 


Ng  i  2  L  /(r  +  A) 

y  ^ 
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with  these  assumptions,  the  Fisher  information  matrix  J  has  the  form 


where 


a  =  L  L  +  3r^  -  ^  (L  +  L  ) 

X  y  2  X  y 

In  order  to  compute  probability  of  anomaly  it  will  be  assumed  that 
the  estimates  of  the  target's  size  and  location  are  constrained  to  be  integer 
multiples  of  the  pixels  or  detector  size,  r. 

Then  the  likelihood  given  by  Equation  (6-17)  reduces  to 


A  = 


EE 

i,j€  target 


(6-28) 


Noting  that 


EE  ''i> 


i,j«  target 


♦ 


where 

=  number  of  pixels  on  target 
one  can  then  reduce  Equation  (6-28)  to 


(6-29) 
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When  the  target  detector  (Equation  (6-29))  is  locked  on  the  target,  X  is  defined 
to  be 


X  =  S{A) 


(6-30) 


which  reduces  to 


X  = 


1  I  AT 

2  N 


o 


2  2 
r  (r  +  A)"^ 

L  L 
X  y 


(nj2  +  NjNg  + 


where  Ng  and  Nj  are  given  by  Equation  (6-24).  For  a  region  of  N.,.  pixels, 
the  probability  that  this  region  will  result  in  an  anomaly  is 


=  P(|  z  I  >  y[Z\) 


=  2  erfc,j,  {yJTX) 


The  probability  of  anomaly  is  then  obtained  by  applying  the  results  from 
reference  1. 


P(anomaly)  <  2 


(FOVj^ 


erfc  (^2X) 
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Plots  of  error  estimation  bounds  and  probability  of  anon.aly  are  given 
in  Figure  6-1  for  the  special  case  of  sampling  at  every  r  unit  (A  =  0).  Also 
included  in  the  plots  are  the  results  obtained  in  reference  1  for  segmentation 
of  a  rectangle  which  has  been  convolved  with  a  blurring  function.  The  close 
agreement  between  the  two  sets  of  plots  for  SNR  2  1.  5  suggest  that  for  A  =  0, 
the  blurring  function  previously  used  can  also  be  used  to  model  the  effects  of 
spatial  quantization. 


7.0  ANALYSIS  OF  ATAC  SEGMENTATION 


7,1  INTRODUCTION 

The  analyses  presented  in  the  preceding  sections  involved  finding 
performance  bounds  for  optimum  object  segmentation.  In  this  section, 
analytic  techniques  are  developed  for  the  specific  segmentor  used  in  the 
ATAC  Autocuer.  There  are  two  main  steps  in  this  segmentor:  edge  detec¬ 
tion  and  interior  point  assignment. 

The  initio.l  step  in  the  ATAC  segmentation  routine  is  the  generation 
of  an  edge  value  map  within  the  window  containing  a  region  of  interest.  For 
this  purpose,  a  Sobel  operator  of  3  x  3  pixels  in  size  is  used.  It  is  assumed 
that  the  (Sobel)  edge  values  corresponding  to  targets  are  among  the  largest 
edge  values  in  the  region  of  interest.  As  used  in  this  algorithm,  an  object 
of  interest  is  a  contiguous  area  consisting  of  interior  points.  An  interior 
point  is  defined  as  a  pixel  that  is  surrounded  in  most  directions  by  edge 
points;  and  an  edge  point  is  defined  as  a  point  with  an  associated  edge  value 
whicli  is  above  some  threshold  value. 

Edge  thresliold  \  alue  is  determined  adaptively  using  the  a\erage  of 
edge  values  within  the  region  of  interest. 

The  1  ule  is 


where 


T  =  K  •  S 


T  =  the  edge  threshold 

S  =  the  average  Sobel  operator  value  in  the  wincfaw 
K  =  a  constant  of  proportionality  (typically  2.Z5) 
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In  this  nnanner  the  extraneous  {not  defining  the  boundary  between 
target  and  background)  edges  are  effectively  discriminated  against,  since 
these  have  the  greatest  contribution  to  the  determination  of  S.  This  edge 
thresholding  technique  is  invariant  to  brightness  and  contrast  variations  in 
the  image  and  has  been  shown  to  give  a  good  target  outline  over  a  fairly 
wide  variety  of  images. 

Once  the  above  threshold  edge  points  have  been  determined,  the 
detection  of  interior  points  within  the  region  of  interest  takes  place.  The 
selection  of  an  interior  point  can  be  easily  visualized.  Because  an  interior 
point  is  surrounded  in  most  directions  by  an  edge  point,  each  ;  ixel  in  tlie 
region  of  interest  is  tested  for  boundedness  by  edge  points  in  eight  direc¬ 
tions  gi\'en  by: 


n360°  _  ,  „ 

0^  “  g  ,n— 1,  .,,,8 


(7-1) 


If  a  pixel  is  bounded  in  6  of  the  8  directions  by  edge  points,  then 
that  pixel  is  labeled  to  be  interior  to  some  boundary.  The  resultant  image 
can  then  be  thresholded  at  6  (for  6  out  of  8  directions)  to  produce  the  interior 
point  iniage.  A  thinning  and  filling  operation  then  removes  isolated  interior 
points  and  fills  in  gaps  in  the  binary  interior  point  image.  The  thinning  and 
filling  operator  consists  of  a  3  x  3  pixel  overlapping  window,  which  is  slid 
within  t!ie  region  of  interest  in  both  the  x  and  y  directions.  If  5  out  of  9 
pixels  in  t!ie  3x3  window  are  occujiied  by  interior  points,  then  the  center 
element  of  this  window  is  set  to  1;  otherw'ise,  it  is  set  to  0. 

Because  of  the  coarseness  of  the  8-direction  edge  point  search,  con¬ 
cave  and  convex  portions  of  objects  tend  to  be  filled  in,  or  shadowed;  tlie 
same  is  true  for  regions  between  two  objects  which  are  in  close  proximity 
to  each  other.  To  overcome  this  drawback,  an  intensity  assignment  of 
interior  cei-sus  exterior  points  is  ca  ried  out. 

To  perform  this  function  the  mean  and  standard  deviation  of  interior 
points  (Up  iTj)  and  exterior  points  o  jr)  computed.  The  assumption 
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is  made  that  the  distribution  of  interior  and  exterior  intensities  are 
approximately  Gaussian.  Tht?  desired  threshold  (between  Uj  and 
is  then  given  by 


T 


<Te  +  aj 


(7-2) 


If  Ug  -  points  >  T  are  selected,  otherwise  points  ■'  T  are  selected.  This 
enables  botii  hot  and  cold  targets  to  be  segmented. 

Following  the  intensity  assignment,  a  determination  is  made  (based 
on  object  area  and  standard  de\'iation)  whether  to  select  the  interior  point 
sketch  or  the  thresholded  in\age  as  the  best  outline  for  the  tai  get. 

Follovsing  this  selection  a  final  thinning  and  filling  operation  is  per¬ 
formed,  to  eliminate  isolated  points  and  to  fill  in  gaps. 

Subject  to  the  thinning  and  filling  operation,  a  connectivity  routine 
is  used  to  identify  distinct  interior  regions  that  are  connected. 

Finally,  a  comparison  of  the  connected  interior  point  regions  thus 
derived  is  made  with  tlie  original  interior  points  in  the  interest  window. 

Those  areas  liaving  the  most  coincidence  of  interior  points  with  edge  points 
are  designated  as  segmented  objects,  and  its  corresponding  original  intensity 
values  are  substituted  into  the  segmented  objects.  Depending  on  the  scenario 
e.xpected  in  the  images,  one  or  more  of  the  segmented  objects  within  the 
interest  window  can  be  tagged  as  likely  objects.  After  this  stage,  a 
calculation  is  made  (based  on  area  and  object  center)  to  see  if  it  is  likely 
that  portions  of  the  desired  target  lie  outside  the  interest  window.  If  so, 
the  window  is  repositioned  a.nd  segmentation  recomputed.  Finally,  the 
segmented  objects  are  passed  on  to  the  high  level  feature  extractor  for 
classification. 

A  statistical  model  for  tliis  operation  and  its  perfoimance  was  devel¬ 
oped  for  a  simple  scene  co  rupted  by  noise.  Then  an  analytic  procedure  for 
finding  the  probability  of  gaps  occuring  within  an  edge  segment  was  developed 
and  verified  experimentally.  The  occurrence  of  gaps  is  important  since  the 
interior  point  assignment  dei^cnds  on  finding  points  surrounded  by  edges.  A 
model  is  then  presented  for  interior  point  assignment,  the  selection  of 
interior  point  threshold,  and  performance  results. 
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7.2  EDGE  ANALYSIS 


As  stated  in  Section  2,  a  simple  model  of  the  imai;e  is  used  to  de\  i'lop 
the  basic  analytic  results.  Since  the  A  TAG  se^;mentor  operates  on  a  wmdov^' 
containing  an  object,  the  model  for  the  scene  is  a  window  N  x  N  pixels  in  si/e 
with  a  background  level  of  zero  average  intensity.  The  object  is  presented  h\ 
an  m  X  m  pixel  square  with  \miform  average  intensity  AT,  centered  in  tie- 
window.  The  entire  window  is  assumed  corrupted  by  Gaussian  noise  ol  zero 
mean  and  variance  ir  . 

The  first  step  is  to  determine  the  distribution  of  edges  as  determined 
by  the  edge  operator,  in  tliis  case  the  Sobel  edge  operator,  defined  by  t:.e 
equations: 

S(x,y)  =  I  Sj^(x,y)j  +js^,(x,y)|  (7-i) 

where 


Sj^(x,  y)  =  I(x  -  1,  y  ,  1)  +  2I(x,  y  -  1)  +  I(x  +  1,  y  -  1) 

-  I(x  -  1,  y  +  1)  -  2I(x,  y  +  1)  -  l(x  +  1,  y  +  1) 


(T-d) 


and 


SY(x,y)  =  Ifx  -  1,  y  -  1)  -f  2I(x  -  1,  y)  -f  I(x  -  1,  y  4  1) 

-  I(x  +  1,  y  -  1)  -  2[(x  T  1,  y|  -  I(x  t  1,  y  '  ll 


(7-5) 


Here  I(x,  y)  represents  the  intensity  of  the  pixel  with  coordinates  x  and  y. 

The  average  values  for  Sj^(x,y)  and  S^(x,y)  are  zero  for  pixels  inside 
or  outside  the  boundary,  but  not  for  pixels  adjacent  to  it. 

tor  pixels  adjacent  to  tlie  boundary  of  the  ohjeet,  eitiier  inside  or  out¬ 
side. 


<^  (Sj^(x,  y))  =  ±4  AT 


(7-h) 
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and 


(S„{x,y))  =  ±4  AT 


(7-7) 


where  denotes  the  ex])ectation  operator. 

If  the  variation  fronithe  mean  is  defined  as 


•^^j|(xiy)  =  -  <'(S,,(x,y)) 


T1 


II' 
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it  is  easy  to  sho\^  tiiiU  the  \  ari;;r>ee  of  Sjj(x,  y)  equals  1  Z  a  .  That  is. 


y) 


.  2 
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similar  result  applies  to  Sy(x,  y).  Furthermore,  ASj_j(x,y)  and  ASy(x,y) 
are  uneo  r  re  lat  ed,  i.e.. 


“^^j|(x,\)  AS^,(x,  y) 


2  2  2  ^2. 

=  a  -a  -a  +o-  =0 


(7-10) 


For  non -liounde  r\  pixels,  the  probability  Pj^.(T)  that  a  pixel  intensity 
doi's  not  exceed  a  yiN'en  tiiresiiold  value  T  is  yi\('n  by 


d 


i\(n  = 


/  ‘'s(x,y 

r\ 


(u)  du 
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'.hero  ti'.L'  pi  ohalnlity  density  function  yj  he  shown  (usiny  NBS  Hand¬ 

book  oi  M  atlie  n  lat  ic  al  f  uiu  tinns  7.  4.  32,  F30i)  to  Ije 


‘^S(x,y) 


> 

4,^ 

o  r  1 

\  2y  ) 

(7-12) 


t  -  ^ 


where  the  error  function  is  defined  as 


and 


2 

cr 


12  <T 


2 


Then 


(7-13) 
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and  the  probability  of  a  detecting  a  false  edge  's  given  by 


- 

.  1 

Pp  " 

/  Ps(x,y)  ^  - 

‘'t 

erf  I 

- 1 

Also,  the  mean  value  for  non-boimdary  pixels  is  given  by 


(7-15) 


^  =  I  ^  ^"S(x,y) 


=  s  =  -iyl^ • 


(7-16) 


Appendix  B  gives  the  analysis  for  boundary  pixels.  It  is  shown  there 
that  t!ie  probability  for  an  tvlge  pixel  being  detected  is  gi\  en  by 


D 


[ 


P„,  .  (a)  du 

S(x,  y) 


(7-17) 


=  1  -  T 


"  (^) 


(7-18) 
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where  the  probability  density  function  for  this  case  is  given  by 


S(x.  y) 


(u)  = 


Zyfrg:  [ 


-{u+4AT)^/4t^ 


+  e 


-(u-4AT)^/4/ 


(7-19) 


Note  Lliat  Equation  (7-18)  for  reduces  to  Equation  (7-15)  for  P^, 
when  AT  =  0.  Similarly,  the  probability  density  function  with  an  edge  present. 
Equation  (7-15),  reduces  to  the  no-edge  distribution  when  AT  =  0, 

Cur^'es  of  the  no-edge  distribution  and  edge  distribution  for  two  signal- 
to-noise  ratios  are  given  in  Figure  7-1.  Also  shown  is  a  displaced  Gaussian 
approximation  to  one  of  the  curves.  It  is  seen  that  an  excellent  approxima¬ 
tion  is  obtained  even  though  the  signal-to-noise  ratio  is  not  extremely  large. 
The  displaced  Gaussian  approximation  is  derived  in  Appendix  C,  where  it  is 
shown  that  the  mean  intensity  for  an  'edge  pixel'  is  given  by  (lor  sufficiently 
large  AT) 


S  =  4  AT 


with  a  variance  gi\  en  by 


(T 


Z 

2 


(7-20) 


(7-21) 


as  befo  re , 

The  second  term  in  Flquation  (7-20)  shows  that  the  presence  of  noise 
causes  the  average  value  to  shift  from  its  value  lAT  in  the  absence  of  noise. 
This  is  a  consequence  of  the  non-linearity  of  the  Sobel  operatoi-. 
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Figure  Distribution  of  background  and  object  edges. 


From  Eq.  (T-lo),  the  moan  \'aluo  of  the  background  and  intorioi- 
pixel  is  4'^^  j  ,  uliilo  from  Eq.  (i-20)  the  mean  value  of  tiie  edqo.s  on  the 
object  boundary  is  4A  1  •*  2 a.  Iho  bobeJ  operator  will  produce  t^^■o  edges 
for  each  true  edge  in  the  picture,  hence  the  average  c'alue  of  tlie  edge-  is 


S  = 


(V-22) 


The  edge  threshold  chosen  is  KS  where  K  =  2.25,  therefore  as  a 
minimum 


r  >  9 


(7-23) 
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so  that  by  evaluating  Eq.  (7-15),  the  probability  that  a  background  pixel  is 
incorrectly  labeled  as  an  edge  is  given  by 

Pp  <  0.022.  (7-24) 

The  probability  that  a  boundary  pixel  is  correctly  labeled  is 
evaluated  using  Eq,  (7-18)  where  the  threshold  is  given  by 

T  =  2.25  S  (7-25) 

Pj^  is  plotted  in  Figure  7-2  for  three  values  of  target  size. 

In  order  to  determine  the  adequacy  of  the  Sobel  operator  in  extracting 
the  square  boundary  in  the  presence  of  noise,  more  needs  to  be  known  than 
tile  probability  that  a  given  edge  pixel  is  detected.  Altliough  P^^  may  be 
high,  tile  probability  ^AD  that  all  edge  pixels  in  the  square  are  detected, 

Pad  ’  ‘Pd'®”  (7-20 

is  smaller  the  larger  the  number  ot  boundary  pixels. 
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Figure  7-2,  Plot  of  edge  detection  probability  vs. 
bNR  for  different  target  sites. 
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Fortunately,  it  is  not  necessary  to  detect  all  edge  pixels  for  adequate 
segmentation  since  a  later  ope  ration  may  still  correctly  separate  interior  and 
exterior  pixels.  Obviously  tlie  least  serious  case  is  wliere  only  a  single  edge 
pixel  goes  undetected,  but  its  neighbors  are  detected.  Examples  of  this  are 
shown  as  A  and  B  in  Figure  7-3.  (To  avoid  unnecessary  computation,  corner 
pixels  such  as  C  and  D  are  excluded  from  the  analysis  since  their  values  in 
the  absence  of  noise  are  6AT  and  2AT,  respectively,  compared  with  the  other 
edge  pixels  whose  values  are  4AT.  Also,  the  correlation  with  their  neighbors 
is  different  for  the  corner  pixels.) 

The  next  most  serious  case  is  where  a  pair  of  neighboring  edge 
pixels  goes  undetected.  Examples  of  this  are  given  in  E,  F,  and  G.  Tlie 
case  of  G  must  be  distinguislied  from  that  of  E  and  F  since  G  causes  a  gap 


Figure  7  3.  Cases  of  interest  for  gap  analysis. 
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in  the  boundary  while  E  and  F  don't.  (Again,  case  N  is  excluded  from  our 
analysis  since  it  involves  a  corner  pixel.)  We  will  call  G  a  first-order  gap. 

Cases  H  and  1  are  second-order  gaps,  while  J,  K,  L,  and  M  are 
third-,  fourth-,  fifth-,  and  sixth-order  gaps,  respectively.  (Note  that  case 
K  is  not  considered  a  fifth-order  gap  since  the  corner  pixel  is  excluded  from 
our  analysis.  Note  also  that  a  gap  of  given  order  also,  in  a  sense,  contains 
gaps  of  lowf  r  order.) 

Before  beginning  the  analysis  three  more  terms  should  be  defined: 
"point",  "gap  point",  and  "non-gap  point".  A  point  is  a  pair  of  neighboring 
pixels  lying  directly  across  the  boundary  from  each  other.  A  "gap  point"  is 
a  point,  both  of  whose  pixels  are  below  the  threshold  for  edge  detection.  A 
non-gap  point  is  a  point  which  is  not  a  gap  point.  An  nth-order  gap  is  then 
formally  defined  as  'n'  adjacent  gap  points. 

A  key  pr.?bability  that  is  needed  is  the  probability  oi  an  ■  h  order 

gap  occurring  in  n  specific  adjacent  points.  By  excluding  corners  fr  gaps, 
this  probability  is  the  same  for  all  possible  n  adjacent  points  in  the  .re. 
For  a  square  of  length  m  there  are  4(m-n-l)  such  configuration  ,  wliere  an 
nth-order  gap  can  occur.  Defining  Dp(n)  as  the  probability  of  detecting  a 
square  not  containing  any  nth  order  gaps  (or  larger),  clearly 

Pj3(n)  =  (l  -  P,^(n))  (7-27) 

The  value  of  n  to  be  chosen  will  in  general  depend  upon  the  siite  of  the  square 
m,  as  well  as  tlic  effectiveness  of  subsequent  operations  in  the  segmentation 
process  in  filling  in  gaps. 

To  obtain  a  formula  for  P^^(n),  define  the  conditional  probability  Pqq 
as  the  probability  that  by  extending  an  nth-order  gap  to  its  neighboring  point 
an  (n  +  1)  -order  gap  occurs. 

+  1)  =  P(^(n)  P^(~  .  (7-28) 
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In  general  Pqq  also  will  depend  upon  n.  However  making  the  assumption 
that  the  correlation  is  small  and  these  can  be  neglected,  will  be  inde¬ 

pendent  of  n.  Experimental  results  discussed  below,  as  well  as  theoretical 
calculations  of  correlation  coefficients,  suggest  that  this  is  a  valid  assump¬ 
tion.  In  Appendix  F,  the  Sobel  operator  correlation  coefficient  between 
next-nearest  neighbors  along  the  boundary  is  found  to  be  only  0,  18  compared 
with  the  value  of  0.49  found  in  Appendix  D  for  nearest  neighbors.  This  high 
value  for  nearest  neighbor  enhances  the  likelihood  of  occurrence  of  high- 
order  gaps  over  what  would  occur  if  there  we  ‘e  no  correlations. 

In  Appendix  E,  the  correlation  coefficient  for  pixels  across  the  boun¬ 
dary  is  found  to  be  only  0.  11.  Consequently,  the  probabiHt'^  of  non-detection 
for  this  case  is  given  to  good  approximation  by  Eqs.  (E27)  and  (E28) 


P  =  P  ‘ 
2S  N 


(1  -  Pd)  . 


(7-29) 


where  is  the  probability  of  not  detecting  any  specified  edge  pixel.  From 
its  definition,  ^23  probability  of  a  first-order  gap. 


Pcd)  =  P2S 


(7-30) 


with  repeated  application  of  (7-28), 


^2S 


(7-31) 


The  conditional  probability  Pqq  obtained  by  specializing  to  the 
case  n  =  2, 


P^(2) 


(7-32) 
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Since  neighboring  pixels  across  the  boundary  are  nearly  uncorrelated, 
it  is  reasonable  to  assume  that  the  same  is  true  for  boundary-parallel  pairs 
of  such  pixels. 


Pc  (2)  =  (P^p)^ 


(7-33) 


where  P2P  probability  of  non-detection  for  a  parallel-to-the-boundary 

pair.  In  Appendix  D,  it  is  shown  how  P2P  calculated  from  certain  graphed 
integrals  of  the  bivariate  normal  probability  function. 

Substituting  (7-33)  and  (7-29)  into  (7-32),  one  obtains 


P 


GC  " 


(7-34) 


Having  calculated  Pq(~;»  rnay  be  calculated  from  (7-28)  or  (7-31) 

Having  calculated  P^(n),  PQ(n)  may  be  calculated  from  P9.  (7-27). 

Experiments  were  conducted  to  verify  the  analytic  techniques 
using  computer  simulation. 

Ten  Monte  Carlo  experiments  were  run  for  each  of  two  target/ 
window  sizes: 

1.  A  16  pixel  square  in  a  24  pixel  window, 

M/N‘^  .2  16/24^  =  1/36; 

2.  A  6  pixel  square  in  a  12  pixel  window 

M/N^  =  6/12^  =  1/24. 

Experimental  averages  and  standard  deeiations  obtained  are  tabulated  in 
in  Table  7-1.  It  may  be  seen  tiiat  the  agreement  with  theory  is  quite  good. 

In  all  cases  the  difference  between  theory  and  experiment  is  well  within  the 
standard  deviation.  The  staiidard  deviations  are  generally  larger  for  the 
M/N^  =  1/24  case  since  fewer  sample  points  wore  involved.  For  tlie  same 
reason,  the  percent  standard  deviation  becomes  larger  for  higher  gap  size. 
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TABLE  7-1.  THEORETICAL  AND  EXPERIMENTAL  VALUES  FOR  SNR  =  4. 
EXPERIMENTAL  VALUES  ARE  RESULTS  OF  10  MONTE  CARLO 
SIMULATIONS  RUN  FOR  EACH  TARGET  SIZE 


M/N^  =  1/36  (=  16/24^) 

M/N^  =  1/24  (=  6/12^) 

Theo  ry 

Experiment 

Theory 

Expe  riment 

s 

84.68 

84. 04  ±  Z.  54 

99.  14 

98,  18  ±  5. 72 

T 

190. 53 

189.0  ±  5.36 

223. 06 

220. 8  ±  12. 58 

^N 

0.  53 

0. 53  ±  0.08 

0.  81 

0. 79  ±  0.06 

^'‘zp 

0.  36 

0,  34  ±  0.07 

0,  70 

0.65  ±  0. 12 

^2S 

0.  Z8 

0. 26  ±  0.06 

0.  66 

0, 59  ±  0. 13 

Pc  (2) 

0.  13 

0.  12  ±  0. 05 

0.49 

0.43  ±  0. 18 

Pa(3) 

0.  06 

0.05  ±  0.03 

0.  37 

0.  33  ±  0.21 

PqU) 

0.  03 

0.02  ±  0.03 

0.  27 

0. 23  ±  0. 22 

Pg(5) 

0.01 

0.01  ±  0.02 

0.  20 

— 

Pc«i) 

0.006 

0.003  ±  0.009 

_ 

0.  15 

— 

Note  that  the  experimental  averages  exhibit  a  bias,  since  all  of  them 
are  less  tlian  or  equal  to  tlie  tlieoretical  values.  The  probable  explanation 
for  this  behavior  is  that,  due  to  cliance,  tlie  average  experimental  intensity 
S  was  (for  both  cases)  smaller  than  the  theoretical  value.  Since  the  threshold 
for  each  experiment  was  set  at  Z.Z5  times  the  average  intensity  for  that 
experiment,  the  average  threshold  was  also  smaller  than  the  theoretical 
v'alue.  (Note,  liowever,  that  due  to  large  roundoff  error,  due  to  use  of 
integer  intensities,  the  experimental  averages  and  standard  deviation  for 
T  do  not  exactly  equal  2.Z5  times  the  cor  respontling  values  for  S.  )  Since  a 
lower  threshold  increases  tlie  probability  of  edge  detection,  and  thus  decreases 
the  probability  of  non-detection,  the  experimental  average  probabilities  will  be 
expected  to  be  below  their  theoretical  values.  It  is  likely  that  even  better 
agreement  with  tlieory  could  be  oiitained  if  the  aposteriori  average  thresholds 
of  Tabic  7-1  were  used  in  making  t.he  tlieoretical  calculations  of  probabilities. 
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Table  7-2  indicates  the  method  of  calculating  the  various  probabilities 
using  the  case  of  the  assumed  edge  detection  of  the  16  x  16  square  in  Fig¬ 
ure  7-3  as  an  illustrative  example.  The  interested  reader  may  check  that 
the  number  of  sample  points  and  number  of  occurrences  is  correct  for  each 
case.  Note  ihat  corner  pixels  have  been  excluded  from  consideration. 

In  Table  7-3  are  tabulated  experimental  and  theoretical  values  of 

P_.(n),  the  probability  of  detecting  a  square  containing  no  nth-order  (or 

^  2 

higher)  gap  for  SNR  =  4  and  object  size  M/N  =  l/36.  The  theoretical  values 

have  been  calculated  from  Fq.  (7-27),  while  the  experimental  values  are 

based  upon  the  10  computer  runs.  The  results  are  satisfactory  in  view  of 

the  small  sample  size.  Again,  the  tendency  of  the  exjje rimental  Pj^(n)  to 

exceed  the  theoretical  \  alue  is  likely  attributable  to  the  bias  in  the  threshold 

which  caused  fewer  gaps  of  a  given  order  to  be  generated,  and  thus  increases 

the  probability  that  a  square  would  have  no  gap  of  this  order. 

Obviously,  it  would  be  desirable  to  make  more  comjjuter  runs  to 
improve  the  statistics,  and  to  make  calculations  for  other  SNRs  of  interest. 

7.3  iriTERIOR  POINTS 

The  second  major  operation  in  the  ATAC  segmentor  is  the  detection 
or  assignment  of  object  interior  points.  Again,  the  simple  scene  model  of 
Figure  7-4  is  used. 

Assume  for  the  object  shown  in  Figure  7-4  tliat  due  to  SNR  reasons 
the  probability  of  an  edge  point  beiiig  correctly  labeled  is  and  the  jjroba- 
bility  for  a  nonedge  point  being  incorrectly  labeled  is  Fp-  Then  the  probability 
that  a  point  interior  to  tlie  object  is  surrounded  by  edges  in  at  least  6  out  of 
8  directions  is  bounded  by 


8Pd(1 


(1  -  1^ 


(7-35) 


A  plot  of  this  function  versus  Pq  is  shown  in  P'igure  7-5. 

For  the  probability'  tliat  a  point  non- interior  to  the  surface  being 
labeled  an  interior  point  consider  the  situations  described  in  Figure  7-(). 
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TABLE  7-2.  TABLE  ILLUSTRATING  CALCULATION  OF  FkO LABILITIES 

BASED  UPON  FIGURE  7-3 


TABi.E  7-3.  EXPERIMENTAL  AND  THEORETICAL  RESULTS  FOR  PD{n), 
THE  PROBABILITY  OF  DETECTING  A  SQUARE  CONTAINING  NO  NTH- 
ORDER  GAP,  FOR  SNR  =  4  AND  M/n2=  1/36  (EXPERIMENTAL  VALUES 
BASED  UPON  10  MONTE  CARLO  SIMULATIONS.) 


No.  oi  Squares 
No 

ntli  -O  rde  r  G  a]} 


PD(n) 


Ex|)o  rimental 


0.2  ±0,4 
0.5  ±0.5 
0. 8  ±  0. 4 
0.9  ±  0.  3 


Theoretical 
10-10 
0. 0007 
0.  04 
0.  2 
0.  5 
0.7 


Figure  7-4.  Simplified  model  for  object  in  noisy  background. 


Case  1  —  Point  on  the  Outer  Boundary  of  Oljject 

In  this  case,  the  point  is  already  surrounded  by  5  directions  as  shown 
in  Figure  7-6.  For  a  6tli  direction  to  be  found  one  point  in  tlie  reniainini; 

3  directions  must  be  found.  This  will  be  bounded  by 


p 

CASE 


1  -  (1  - 


(7-36) 


Using  N  =  25  and  P^.  =  0.02  from  (7-0)  this  yields  approximately  51  percent. 


Case  2  —  Point  Outside  the  Outer  Ilottndary  but  Closer  than  m  i^ixels  to  the 
Object 

lli'i'e  thi’  pixil  is  surrounded  in  three  directions,  to  be  labeled  as  an 
interior  point  it  must  be  surrounded  by  at  least  three  in  the  remaining  five 
directions  —  this  is  bounded  by 


CASE  2 


(1  -  -t- 


(1  -  q)  -I- 


(7-37) 


f 


T}ic  thinning  operation,  which  requires  5/9  adjacent  cells  to  be 
interior  points  prior  to  assigning  an  interior  point  will  therefore  remo'/e 
virtually  all  of  the  misassigned  points  in  Cases  3  and  2,  and  will  remove 
tile  majority  of  the  points  assigned  in  Case  1. 


Intensity  Assignment 

from  the  preceding  discussion,  tiie  measured  interior  jioint  average 
intensity  will  be’I' 


AT  X  m'^  -t  0  X  (4m  t  4) 
(m  +  2)“ 


(7-40) 


INTAV 


=  AT 


\m  +  2/ 


The  measured  exterior  point  a'.erago  will  be 


LiiXT.lV 


.  =  0 


{7-.-, 


(7-42) 


Ll,. 


The  measured  exterior  point  standard  dexiation  will  be  fairly  close 


ilX  f 


=  cr 


(7-43) 


And  the  mensui'ed  interior  point  stauflarrl  deviation  will  lie 


(7-44) 


This  assumes  that  alt  of  ti.e  boundarv  points  ;ire  Ctise  ’  points  and 
not  reiuoxwd  by  thinniu;;. 


The  selected  th.reshold  will  be 


(7-45) 


(7-46) 


is  a  function  of  object  si/e.  In  practice  when  m  <  8,  (resulting  in  small  p  ), 
the  interior  points  pictui-e  is  selected.  As  p  — '  1  this  threshold  approaches 
the  optimal  thresliold  of  AT/2.  As  the  SNR  (AT/ct)  is  increased,  the  chosen 
thresliold  becomes  increasingly  tolerant  (lower),  allowing  more  weight  to  the 
inclusion  of  oliject  points  than  to  tlie  exclusion  of  unwanted  noise  points, 
since  tliese  will  be  cleaned  up  in  later  thinning  and  connectivity  stages. 

Following  the  selection  of  tlu'  threshold,  the  window  is  passed  through 
tliis  threshold.  D/pending  on  SNR  this  tlireshold  will  determine  tlie  Nqq  = 
Probability  of  correct  assignment  of  object  points  and  =  probability 

of  correct  assignment  of  background  points.  A  plot  of  these  probabilities 
is  shown  in  Figure  7-7.  For  a  connected  object  to  maintain  its  ijroper  sliaj^ie 
to  within  5  percent,  Pqq  ^  0.95  and  P^^j^  >  =  1  -  l/2  m  (for  few  u;i\'- anted 
points  to  exist  on  the  boundary  of  the  extracted  object)  tlu;  carves  of 
Figttre  7-7  show  an  SNR  requirement  of  about  5  for  T  =  A  F/ 2  and  m  =  10. 
This  requirement  is  consistent  with  the  requirements  on  tlie  otiier  stages 
of  this  process. 
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7.4  SUMMARY  AND  CONCLUSIONS 


Of  course,  if  the  only  problem  were  to  extract  the  shapes  of  simple 
objects  in  idealized  backgrounds  such  as  that  illustrated  in  Figure  7-4, 
there  would  be  a  number  of  simpler  techniques  which  would  be  capable  of 
performing  an  effective  segmentation.  The  techniques  described  here,  are, 
however,  applicable  over  a  wide  range  of  circumstances  encountered  in 
actual  IR  images  and  not  described  by  the  model  of  Figure  7-4.  It  is  the 
purpose  of  this  analysis,  however,  to  show  in  broad  terms  the  noise 
sensitivity  of  the  algorithms  that  are  employed. 
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8.  0  CONCLUSIONS  AND  RECOMMENDATIONS 


An  analytic  modeling  approach  for  image  processing  algorithms  has 
been  developed  and  applied  to  the  complex  object  segmentation  operation  used 
for  automatic  target  classification/recognition  in  FLIR  images. 

Starting  with  a  simple  scene  model,  performance  bounds  were 
defined  for  the  segmentation  operation  and  analytic  means  for  evaluating 
these  bounds  were  developed,  A  technique  was  developed  to  extend  these 
results  to  more  complex  scenes  involvung  more  detailed  target  shapes  and 
presence  of  simple  clutter,  while  maintaining  the  tractability  of  the  analytic 
evaluation  of  performance  bounds.  These  results  will  provitie  the  means 
for  obtaining  useful  analytic  results  for  several  problems.  First  the  per¬ 
formance  of  a  specific  segmentation  algorithm  on  images  with  noise  and 
simple  clutter  can  be  compared  with  the  performance  bounds  to  determine 
if  that  segmcntor  is  performing  as  well  as  can  be  expected  or  is  deficient 
in  some  respect,  e.  g, ,  in  responding  to  noise  or  the  presence  of  clutter. 

The  analytic  techniques  can  be  used  to  extend  the  complexity  of  target 
and  clutter  description  to  find  the  degree  of  complexity  necessary  for  good 
analytic  results  and  the  important  descriptive  features  for  the  particular 
image  processing  operation  considered.  This  approach  can  be  used,  e.g., 
to  find  a  relatively  simple  analytic  model  of  a  particular  FLIR  image  that 
can  predict  segmcntor  performance  analytically,  without  the  statistical 
uncertainties  of  testing  a  segmcntor  on  a  small  number  of  images. 

The  use  of  these  performance  bounds  also  will  allow  the  specification 
of  maximum  permissible  system  noise  characteristics  for  achieving  a  spe¬ 
cified  segmentation  performance. 
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In  continuation  ot  tlus  effort  it  is  recommended  that  the  previous 
results  be  extended  with  the  j^oal  of  establishing  and  verifying  analytic 
models  for  tlie  performance  of  target  recognition  and  tracking  systems. 

Tile  recommendations  are  di\ided  into  5  primary  tasks.  Task  1, 
denoted  '  Image  Modeling''  will  extend  image  models  to  include  3-dimension 
effects  and  impro\'e  the  clutter  models.  Task  11,  denoted  "Tracker  Modeling" 
will  extend  the  m  idels  to  include  Target  Tracking  Systems.  Task  III, 
denoted  "Performance  Analysis"  will  involve  derivation  of  improved  perform¬ 
ance  bounds  for  target  recognizers  and  trackers  and  extension  of  these 
bormds  to  include  clutter  limited  performance.  Task  I'V,  denoted  "Context 
Modeling"  will  dec'olop  methods  for  use  and  analysis  of  non-silhouetted 
information,  e.g.,  wlieels,  treads,  engines.  Task  V,  denoted  "Image 
System  Modeling"  will  apply  the  modeling  philosophy  and  the  results  of  the 
pre\'ious  tasks  to  specifically  model  a  target  recognizer  anc’  a  target  tracker. 

Specifically,  in  Task  I; 

1.  Derive  complete  3-D  target  models 

3.  Derive  impro\'ed  3-D  clutter  models 

3.  Model  shape  recognition  for  2-D  pi-ojections  of  3-D  objects. 

4.  Generate  synthetic  2-D  scene  models  for  use  in  performance 
ajialy  .sis. 

In  Task  11: 

1.  Characterize  Target  Trackers 

2.  Develo])  analytic  models  for  Target  Trackers 

3.  Vei'ify  (he  models  using  simulation 

In  111: 

1.  Deri\'e  imj^roveu  pe  ifo  T-mance  l^ounds  for  target  recognizers 

2.  Derive  performance  boimds  for  Target  Trackers 

3.  Tlxtend  our  performant:e  bound  apjjroacli  to  inchide  clutter 
limited  (as  opposed  to  noise  limited)  effects 

4.  Kstablisli  the  limiting  factors  on  target  tracker  performance 

5.  Develop  ajialyticnl  meajis  of  pi-edicting  performance  of  trackers 
eind  recogni  ze  r  s 
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In  Task  IV: 


1.  Derive  models  for  typical  target  contextural  parts:  wheels, 
treads,  engines 

2.  Establish  a  model  for  use  of  contextual  information,  both 
individually  and  synergistically 

3.  Develop  performance  models  for  the  use  of  this  information  in 
target  recognition. 

In  Task  V: 

1.  Implement  a  target  recognizer  model  using  the  image  meta-model 
philosophy 

2.  Implement  a  target  tracker  with  the  meta-model 

3.  Verify  the  accuracy  of  the  models  implemented  in  (1)  and  (2), 
above . 
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APPENDIX  A 


New  approach  to  FLIR  segmentation 

lawrence  M.  Rubin  and  Richard  L.  Frey 

Electronics  Department,  Hughes  Aircraft  Company 
Culver  City 


Abstract 


An  approach  to  the  analysis  of  shape  extraction  (segmentation)  operators  is  presented 
based  on  maximum  likelihood  parameter  estimation.  A  simple  image  model  is  presented  and 
objective  measures  of  segmentation  quality  are  presented.  Bounds  on  the  noise  performance 
of  any  segmentation  operator  are  derived  and  these  are  compared  with  experimental  results 
obtained  on  sample  images. 


Introduction 


In  systems  requiring  pattern  recognition  techniques  to  be  applied  to  imagery  such  as 
that  generated  by  Forward  Looking  Infra-Red  (FLIR)  sensors,  a  critical  stage  is  often  the 
extraction  of  object  shapes  from  the  scene.  These  shapes  can  then  be  measured  and  the 
resulting  features  applied  to  a  classifier  for  such  processes  as  automatic  recognition  of 
gargets.  A  variety  of  algorithms  have  been  developed  to  perform  this  shape  extraction 
function'"'' .  These  algorithms  usually  employ  such  techniques  as  edge  detection,  intensity 
thresholding,  and  some  type  of  connectivity  criteria  to  obtain  their  estimate  of  object 
shape.  In  general  these  algorithms  are  heuristic  in  nature  and  typically  their  performance 
is  judged  subjectively.  Factors  such  as  noise  sensitivity  are  usually  not  brought  into 
play  in  the  evaluation  of  these  techniques.  Useful  analytic  results  pertaining  to  these 
algorithms  are  difficult  to  obtain,  due  to  both  the  difficulty  in  modelling  realistic 
scenes  and  the  intractability  of  the  algorithms  themselves.  This  paper  will  attempt  t.‘ 
study  segmentation  algorithms  on  a  slightly  more  formal  basis,  using  results  obtained  from 
maximum  likelihood  parameter  estimation  theory.  A  simple  image  model  will  be  utilized, 
and  the  effectiveness  of  a  segmentation  procedure  will  be  gauged  in  terms  of  its  accuracy 
in  estimating  the  parameters  of  the  model.  Bounds  will  be  derived  on  the  performance  of 
an  optimal  procedure  and  these  bounds  will  be  compared  to  the  performance  obtained  from 
several  simulations  of  other  segmentation  procedures. 

Image  model 

The  model  that  we  will  employ  to  represent  sample  images  is  illustrated  in  Figure  1. 

We  use  continuous  coordinates  in  both  x  and  y  directions. 

g(x,  y)  is  used  to  represent  the  intensity  distribution  of  the  image 

g!x,y)^b(x,y)*s(x,y)+n(x,y)  (1) 

0  ^  X ,  y  ^  R 


where  b(x,  y)  represents  the  background  intensity  b(x,  y)  -  b. 


s(  X  ,  y  )  =  AT  •$,•(  X  )  •  S 

hy) 

(2) 

S , ( x)  =  rec  t ( X  -  C^, 

) ©  nec t( X ,  r) 

(3) 

and 

S,(y)  =  rect(y  -  , 

Ly) ©  rect(y ,  r) 

(4) 

where Q  i  nd  i  cates  convolution 

1  0  hi 

>  1/2 

(5) 

rect(x,  t)  = 

'l/i|>'l 

1  1/2 

( C  _  C  )  1  s  the  uir’.i-  of  the  object  rectangle  and  L  ,  L  is  the  length  of  the 
rec  t  a  d'l  I  ^  i  n  ttiu  «  and  y  directions  respectively.  r  represents  a  spatial  band  limiting 
process  dti  lied  to  to-  ifage  (see  figure  ?).  AT  represents  the  intensity  amplitude  of  the 
rectangle  and  nix,  y)  represents  a  stationary,  zero  mean  Gaussian  spatially  white  process 
such  that 
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(6) 


E[n(x,  y),  n(x  +  x,,  y  +  y,)]  =  ^  x  Uq{x,,  y.) 
where  Uq( • , • )  is  a  two  dimensional  impulse  function. 

Per f ormance  bounds 

We  state  the  segmentation  problem  as  follows:  Given  an  observation  of  g(x,  y), estimate 
the  center,  size,  intensity  and  shape  of  the  rectangle.  We  will  treat  the  parameters 
C  ,  C  ,  L^,  L  ,  and  aT  as  non-random;  that  is  to  say  with  no  prior  known  underlying  prob¬ 
ability  distribution. 


Under  these  circumstances  the  problem  reduces  to  a  mu  1 1 i - pa  name t er  estimation  problem 
and  results  relating  to  the  bounds  on  performance  of  such  an  estimate  are  well  known'. 

These  results  will  be  repeated  below: 

Let 

g(x,y)=s(x,y,  F)+n(x,y)  (7) 

where  F  is  a  vector  of  parameters  =  [fi.  fj.-.-f^]  that  will  be  referred  to  as  the 
feature  vector  and  n(x,  y)  is  as  defined  in  (6).  The  f-  is  considered  to  be  non-random. 

Let  fj  be  an  unbiased  estimate  of  f^-.  Then 

variance  [f-  -  f^]  >  E^-  (8) 

where  E.-  is  the  t^^  diagonal  element  of  the  matrix  E.  given  by  equations  (9)  and  (10). 

E  =  J-‘  (9) 

where  J  corresponding  to  the  Fisher  information  matrix  is  given  by 

R 


where  R  is  the  region  of  obsor'vation  of  u(x,  y) 


The  interpretation  of  equations  (8)  -  (10)  is  as  follows:  Given  a  scene  containing 
an  object  or  objects  characL erizable  by  a  vector  of  features  F,  any  segmentation  procedure 
intended  to  extract  one  or  "ore  of  the  intended  objects  will  produce  a  shape  from  which  the 
error  variance  cf  any  unbi.iied  estimate  of  the  intended  shape  features  is  bounded  by 
equation  (8). 


lo  show  hOw  Itie  bound  applied  to  the  model  of  section  (II)  we  note  that  letting 
F  M  P I  .  C  ,  C  ,  L  ,  I  ]  ^ 

‘  X  ’  y  X  y 

q(x,  y)  as  given  in  equatums  (1)  -  ( f )  .  and  r-.vL  r<<L 
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After  inverting  J  and  u'.inn  (8),  the  error  bounds  can  be  seen  to  be: 


var[Al  -  AT] 


No 


2  L  L 


var[ 
va  r[  t 


Cl 


c  ] 
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N  O 
4 


U  Cl 

4 


x'-y  1-2  r'/L^L^ 

(r/l,^) 

A  I  ■■ ' 

(  r  /  1  ) 

A 1  ■ 


A-. 


(11) 


(12) 


(13a) 


(13b) 

(13c) 


''ar[Ly  -  L„]  > 


yJ  i  ZTFI 


[4  -  2  r/Lj 

2AT^Ly 

2  - 

r(l/L^  ♦  1/Ly)  - 

*rr- 

-No  n  J 

r 

[4  -  r/L^] 

X  y 

r(l/L^  ^  VL^)  -  4^ 


X  y  J 


Anoma lies 


(13d) 

(13e) 


The  bound  of  equation  (8)  is  derived  from  the  behavior  of  an  estimator  In  the  neighbor¬ 
hood  of  the  niaximuni  likelihood  estimate.  A  further  consideration  is  the  probability  of 
anomaly,  i.e.,  the  probability  that  the  estimate  lies  nowhere  near  the  maximum  likelihood 
estimate. 


The  case  tiiat  will  be  discussed  is  the  probability  of  anomaly  for  maximum  likelihood 
extraction  of  squares.  To  perform  the  analysis  a  sufficient  statistic  for  the  likelihood 
ratio  hill  first  be  derived.  A  bound  on  the  probability  of  a  given  noise  area  having  a 
greater  likelihood  than  the  actual  square  will  be  given,  and  this  bound  will  be  used  to 
bound  the  overall  probability  of  anomaly. 

Following  equations  M)  -  (5),  we  will  set  L  =  L  =  L  (for  a  square),  b  =  0  (this 
assumes  known  background)  and  r'-  L.  To  good  a ppr ox  1  me t i on  then  a  suitable  likelihood 
func  t  ion  to  use  is 


wheie  d  =  L/2 


Letting  G(C^  ,  ,  L) 


y)  dx  dy 


(14) 


(IG) 


the  value  of  .'.1  that  maximizes  (14)  for  any  choice  of  and  L  Is  2(0^^,  C^,  L). 

Under  these  circumstances  (14)  reduces  to 


A.(C  .  C  .  L)  =  A„(C  .  C  .  L  AT)  =  GMC_  C„.  L)  •  L' 

Ay  *  j  I  *■  y 

For  a  square  centered  at  u,v  of  side  I  and  offset  imbedded  in  noise  as  defined 

equation  (6),  the  expected  value  of  this  statistic  taken  about  the  center  of  the  square 

X  =  E[gMu.  V.  =  AT^’l^  e 


(16) 

n 

1  s 

(17) 


where  the  second  term  on  the  right  of  (17)  Is  due  to  noise. 


The  probability  of  G’(C  ,  C  ,  L)'L'  exceeding  X 
L/2,  C,  -  L/2  y  <  C„  ♦  L/?  contains  only  noise  of 
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where  the  region  - 
spectral  amplitude 


L/2  < 
1  s 


(18) 

(18a) 


=  2  erfc' 


Vno' 
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(  IRb) 


whei  e 


erfc* ( x)  i 


_L 

/H 


d2 


(19) 


Note  in  particular  that  the  right  side  of  (18b)  is  independent  of  L.  Therefore,  the 
probability  of  any  square  noise  region  having  a  statistic  greater  than  X  is  independent  of 
the  size  of  the  region.  An  image  of  area  R'  can  be  viewed  as  containing  R'/L^  statistically 
independent  square  regions  of  area  L^.  Therefore,  an  approximation  to  the  probability  of 
an  anomaly  of  size  can  be  computed  using  the  union  bound  as 

p(Anomaly  of  Size  L’)  =  ^  x  2  ®'‘fc*(^^)  (20) 


To  estimate  the  total  probability  of  anomaly  these  probabilities  must  be  summed  over 
the  different  sizes  of  squares.  Before  performing  the  summation  however,  the  correlation 
of  different  size  squares  centered  at  the  same  position  in  the  image  must  be  taken  into 
account.  What  follows  relates  to  th’s  correlation. 


Lemma  1  -  Let  x,  and  Xj  be  zero  mean  statistically  independent  Gaussian  variables  with 
standard  deviations  oj  and  o.:  respectively.  Given  xj  <  Tj  and  Tj  <  T^,  the  probability 
that  X,  +  xj  >  Tj  is  bounded  by 


P(xi tx^vT,  )  < 


■  -.--iiz-  iexp 
Voi  '+0: ' 


T  Z  Z 

.i[L_{i.°j_  )] 
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erfc’ 


^l-erfc*(T,/o, )  j 


Proo  f : 


J _  g'X I ■ /2o  , ' 


P(x,+x,.'Tj  |x,  <1,)-J 


/Ti  n  I 


^1  -erf c* (T , /o ,  ) j 


erfc*  ^(  Tj  -X,  )/oj^  dx, 


(21) 


(22) 


Now  using  the  bound 
erfc*(x)  <  ^e’’' 

for  X  >  0,  we  get  by  substituting  (2d)  into  (22) 


P(x,  +  x,-  Tjlx,-  T,) 


-x,V2ii"  -(Tj-x.)  /2oj' 


/2-  n, 


i>v 


dxi 


^1  -erfc‘{l  ,/-3,  )j 
Completing  the  square  in  the  exponent  of  (24)  obtains 


,  (1-^)] 

P(x,  +  x,>Tj  )  <  - X  Sy - - - 

r^"  n,  ll-erfc*(T,/(i.l  I 
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dx  1 
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(24) 


(25) 


with 


eq 


(25a 
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(25b) 


dnd  M  -  T2  — ^*^7 


From  (25)  then 


o 

P(x,+X2>T2)  <  - - 

/ZtT  O  1 


X  4exp 
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(erfc^y^)) 

^-erfc*{T./o,)^ 


Substituting  (25a)  and  (25b)  into  (26)  gives  (21).  The  above  lemma  can  be  used  in  forming 
the  theorum  below. 

Theorum  1:  Given  the  region 

-  L2/2  <  X  <  Cjj  +  L2/2 

Cy  -  L2/2  <  y  <  Cy  +  L2/2 

to  contain  white  Gaussian  noise  of  spectral  amplitude  No  and  given  Li  <  L2. 


/  \  -A/No  I  2  I  2 

P(L,-G(C^.Cy.L2  )  >  /X  I  L,-G{C^.C  .L,)  <  /\  j  <  he  •  F(  ■  2A/No) 


1 -erfc* I s 


1])) 

[ 1 -erfc*{s )] 


Proof : 


Let  X,  =  G(C^.  C  ,  L,M  ^ 


G(C^.  C  ,  LjXj'  -  G(C^.  C  .  L,)  L,' 


Then  x,  and  Xj  are  tero  mean  Gaussian  random  variables.  Furthermore,  x,  and  X2  can  be 
shown  to  be  statistically  independent  with  standard  deviations 


,|No  U 

.Iho  Vl  2’-L  1 

•  Vt 


res  pec  t  i  ve ly 


Therefore, 


p(l,-G(C^,C^,L  J  >  A  1  L,-G(C^,Cy,L.,  ,  A^=  ,(x,+X2  >  A  1  x .  <  A  ^)  (3?) 

which  from  substitution  into  Lemma  1  and  after  some  simplification  of  terms  yields: 

X  ,e  X  (33) 

*  [ 1 -erf c* ( A /No/2 ) ] 


Making  the  substitution 


s  -  X/(No/2) 


yields 


ri-erfc*fs-  -  -J:  (1- 

_ \  \  1  ♦  W  /W  '  / 1  A  w  ' 


[  1  -  er f c* ( s ) ] 


or 


1  -  erf  c ' 


—  ( 1 


z  [/T+w  -  1]) 


[1 -erf c*( s ) ] 

which  is  the  desired  expression. 


(35) 


Note  that  F(0,  x)  =  0  and  that  F(w,  s)  -►  1  as  w  F(w,  s)  is  a  monotonically 

increasing  function  of  w  representing  the  percentage  of  asymptotic  probability  as  L?  -► 

A  plot  of  F(w,  s)  for  various  values  of  s  is  shown  in  Figure  2.  From  this  figure  it  can 
be  seen  that  for  w  =  1  the  size  probabilities  are  relatively  uncorrelated,  achieving 
approximately  .707  of  their  full  scale  value  for  s(which  represents  total  signal  to  noise 
ratio)  exceeding  about  2.  Therefore,  -  Li^/L,^  =  1,  or  =  2-Li",  squares  of  size 

Lj^can  be  treated  to  have  anomaly  probabilities  approximately  independent  of  size  Lj^ 
squares  on  the  same  centers.  By  doubling  the  area  of  the  previous  sample,  then,  independent 
samples  can  be  obtained. 


Allowing  the 
P(Anonia1y)  < 


minimum  size  square  to  be  r^  and  using  equation  (20)  and  the  union  bound: 
R 


.7  X  2  e  r  f  c  *  ( 


^  ^  X  2  erfc*  {M)  +  ^  x  erfc*  {M) 


(36) 


=  2  erfc*  ^  (1  .  t  .  i  +  ...) 

=  <\£' 


Experimental  resul ts 


To  illustrate  the  bounds  obtained  in  the  previous  section,  tv;o  experimental  image 
operators  were  applied  to  a  series  of  images.  Each  image  measured  100  x  100  pixels^  in 
area  and  each  contained  a  5  x  5  pixel'  square  at  its  center.  The  images  were  corrupted 
by  uncorrelated  noise  of  a  varying  amplitude.  Fifteen  independent  images  were  obtained  at 
each  noise  level  tested.  The  function  of  each  operator  was  to  estimate  the  position,  size 
and  intensity  offset  of  the  square  in  each  image.  A  sample  result  of  such  an  operator  is 
shonn  in  figure  3.  The  errors  in  estimated  parameters  and  probability  of  anomaly  were 
calculated  at  the  varying  noise  levels  and  compared  to  the  computed  bounds  obtained  in  the 
prev  i  ous  sec  t  i  oti . 

The  first  operator,  denoted  OPl,  is  basically  an  exhaustive  search  for  the  maximum 
likelihood  square  over  the  ir.age.  Each  possible  size  square  is  tested  in  every  possible 
position,  and  a  likelihood  meaturi'  similar  to  equation  (16)  is  computed  for  each  combina¬ 
tion.  The  conibination  having  the  greatest  value  of  this  measure  is  then  used  to  estimate 
the  relevant  parameters.  Tiie  second  operator,  denoted  0P2  ,  operates  in  a  similar  manner, 
but  uses  an  iterative  3x3  smoothing  filter  to  estimate  the  average  intensity  of  square 
whose  effective  area  is  a  linear  function  of  the  number  of  smoothing  iterations.  The 
intention  was  to  determine  the  pe r f o mia nc e  degradation  caused  by  replacing  a  global 
averaging  process  by  a  series  of  local  averaging  processes.  The  results  obtained  from 
these  operators  are  graphed  in  Figure  A  which  shows  in  general  a  good  agreement  with  the 
bounds  obtained  in  the  previous  sections. 

Summa  ry  an^  cone  1  us  ions 

In  general,  the  determination  of  noise  sensitivity  of  a  segmentation  operator  over  a 
wide  variety  of  imagery  types  is  a  difficult  procedure.  By  simplifying  the  imaqe  mode, 
however,  bounds  can  be  obtained  for  ttie  error  variances  of  certain  shape  parameters  of 
which  it  IS  often  the  object  of  the  seomentation  procedure  to  estimate.  Real  image 
operators  can  be  compared  with  these  bounds  to  provide  objective  measures  of  performance, 
determine  SNR  deq r a d  .i  1 1  on  ,  and  suggest  potential  areas  for  improvement.  As  important  as 
tne  bounds  on  parameter  cstimaiton  are  the  bounds  on  anomalous  behavior,  which  can  be  used 
to  bound  the  tradeoff  between  detection  probability  and  false  alarm  rate  for  systems 
applying  segmentation  to  functions  such  as  automatic  target  acquisition.  A  simple  image 
model  was  employed  in  this  paper,  and  llie  derived  bounds  were  compared  with  actual 
performance  otita'ned  fro'  two  operators.  The  observed  data  agreed  well  with  the  computed 
bounds.  In  tl.e  future,  we  plan  on  comparing  the  bounds  with  performance  obtained  from 
more  realistic  s  t  g  n'en  t  a  1 1  on  opierators,  as  well  as  extending  the  range  of  a  pp  1  i  c  a  b  i  1  i  ty  of 
the  image  model  er. ployed. 

Z\-C 
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Figure  1 .  Image  Model 
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APPENDIX  B.  SOBEL  OPERATOR  PROBABILITY  DENSITY 


This  appendix  will  derive  the  probability  density  of  the  Sobel  edge 
operator  defined  in  Eqns.  (7-3)  -  (7-3).  To  compete  this  result  it  should 
be  noted  that  the  horizontal  and  vertical  contributions  to  the  Sobel  opt;r.i';j'- 
Sj^(x,y)  and  S^,(x,y)  in  (7-1)  are  uncorrelated  random  \  oriabl  js. 

The  probability  density  function  P(S)  for  the  sum  of  the  absolute 
values  of  two  Gaussian  independent  variabJes  with  non-zero  means  is  given 
by 


P(S) 


dy  6(S  -  X 


y  )Pj.  (y), 

o  ’  o 


(B-1) 


whe  re 


-(x-x  )^/2cr‘^ 
o 


(B-3) 


j  -(v-y^)^/Zg^ 

P  (y>  =  e 


(B-3) 


Zn  (T 


It  is  reaciily  verified  tlial  P(S)  is  properly  normalized 


l'(s)  ds  =  i 


(B-4) 


by  interchanging  the  order  integration. 


B-1 


Breaking  the  integration  in  (B-1)  into  its  contribution  from  the 
four  quadrants  of  the  x-y  plane. 


P(S)  =  /  dx  /  dy  6(s  -  X  -  y)  (x)  (y) 


+  /  dx 


CO 

/  dy  6(s  +  X  -  y)  p  (x)  p  (y) 
/  '  o 

r\ 


r 


dy  6(s  -  X  +  y)  p^  (x)  p  (y) 
o  ^  o 


(B-5) 


dx  /  dy  6(s  +  X  +  y)  p^  (x)  p  (y) 


o  '  o 


Changing  integration  variables  from  x  to  -x  and  y  to  -y  as  appropriate, 
and  using  the  fact  that 


o  o 


(B-6) 


(B-7) 


Eq.  (B-  i)  becomes 


1>(S)  =  /  dx  /  dy  6(s  -  X  -  y)  (x)  +  p_^  (x)  p  (y)  +  p  (y) 

II  0  0^0^0 


"^X  y  +  P-x  ,y  +Px  ,-y  +  P-x  ,  -yj 

0^0  u’ ^o  o’  'o  o  'o 


(B-8) 


B-Z 


where 


P 


X  y 
o  o 


(S) 


d 


dy  6(s  -  X  -  y) 


(B-9) 


Substituting  (B-Z)  and  (B-3.)  into  (B-9),  results  in 


P  (S)  = 

X  y  ^ 

O''  o  Zncr 


CO  < 

b  - 


d-y  6(s  -  X  -  y)  e 


(x-x^)^(y-y^)2 


Changing  to  integration  variables  U  and  V  defined  by 


/2  2"^ 


(B-ll) 


U  =  X  +  y 


(B-12a) 


V  =  X  -  y 


(B-12b) 


o  r 


X  = 


U  +  V 


(B-13a) 


Y  = 


U  -  V 


(B-13b) 


The  Jacobian  of  tlic  transformation  is  equal  to  l/2,  and  since 

2  2 

(X  -  +  (y  -  yf  X  [I  (U  +  V)  -  4  [i  (U  -  V)  - 


!l(U^  +  V^)  +  x^  +  y^  -  U(x  +  y  )  +  v(y  -  x  )| 
(2  o  ^  o  o  '  o  '  o  o  I 


(B-14) 


B-3 


^  / 

—  ^  r\ 


du  (s  -  U) 


dV  exp  '  -(B14)/Z 


(B-15) 


where  B14  is  the  last  expression  in  Equation  B-14. 

Using  the  NBS  Handbook  integral  7.4.32  yields  the  result 


From  (B-8)  it  then  follows  that  the  density  distribution  function 
P(S)  is  given  by 


P(S)  = 


r  2^  2^_2‘1  2 

-  X  +y  +S  /4(r 
1  o  o  J 


4-v/TTcr 


F(S,x  ,y  )  +  F(S,  -X  ,  -y  )  +  F(S,x  -y  ) 
I  '  0  -^0  o  o  o  ’o 


+  F(S,  -X  ,  -y  ) 

'  ’  o  'o  I 


(B-19) 


which  is  the  general  result. 

Specializing  to  the  case  of  a  horizontal  edge  pixel 


X  = 

o 


SyCx.y)  =  0 


{B-20a) 


iS'|^SH(x,y)]  = 


y^  =  <S'|  ST,(x.y)  I  =  4AT 


{B-20b) 


Eq.  (B-19)  reduce.s  to 


P(S) 


(4AT)^+S^ 


/4  (T 


[f(S,  o. 


4AT)  +  F(S,  O,  -4AT) 


(B-21) 


and  (B-18)  reduces  to 


F(S,0,  4AT)  =  e 


_4SAT/22 


u  rf 


fS  +  4AT\ 


2  (T 


+  c  rf 


(B-22) 


B-5 


Thus  tlie  probability  density  of  a  Sobel  edge  operation  in  the  vicinity  of  an 
intensity  discontinuity  of  magnitude  AT  is 


P(S)  = 


^-(S-4AT)^/4£^ 


+ 


^-(S+4AT)^/42^ 


•  erf 


Is  +  4AT 
\  2(t 


+ 


erf 


(B-23) 


which,  aside  from  a  change  of  variable,  is  the  same  as  Eq .  (7-19)  in  the 
main  text. 

From  the  definition  of  the  error  function,  Eq.  (B-2  3)  can  also  be 

written 

2 

(B.24) 


*"(3)  “  I  h 


erf 


/S  H-  4 
\  2(r 


4AT  \ 


+  erf 


(S  -  4AT\ 

^  2^  / 


B-e 


APPENDIX  C.  AVERAGE  VALUE  AND  VARIANCE  OF 
SOBEL  EDGE  OPERATOR 


This  appendix  will  calculate  the  mean  and  variance  of  a  Sobel  edge- 
operator  in  the  vicinity  of  an  intensity  discontinuity  in  an  image.  The 
definition  of  the  Sobel  operator  may  be  rewritten 


S 


Sr  I  +  I^V 


(C-1) 


where  and  Sy  are  the  horizontal  and  vertical  components  respectively. 
The  variance  of  S,  a- is  defined  by  (cf  Eqn.  7-3  -  7-5) 


(^2^=  -  (.g’CS))  ^  . 


(C-2) 


Assume  we  are  dealing  with  a  vertical  edge  with  sufficiently  high 
SNR  that 


f  3y  ^  f  (Sy)  =  4AT 


(C-3) 


Now  Sj_j  can  be  written  as 


<f-  S 


II 


1  1  1 
-  /  ^  ® 


dx 


TT  cr  J 


z 

-  dx 


(C-4) 


C-1 


Hence 

=  {4AT)^  +  2a^  +  8AT\/f  (C-IO) 

Substitvitinj'  (C-10)  and  (C-5)  into  (C-Z)  gives 

(T ^  =  2  £.  »  (C-11) 

v'.iiicli  derives  Eq.  (7-2  1)  of  the  main  text.  Note  that  although  (C-5)  and 
(C-11)  wore  derived  for  the  case  of  a  vertical  edge,  a  horizontal  edge  would 
gj\e  tiio  same  resvilts. 


C-2 


APPENDIX  D.  PROBABILITY  OF  NON-DETECTION  FOR  A 
PARALLEL-TO-THE  BOUNDARY  PAIR 


In  this  appendix  we  calculate  the  probability  P^p  due  to  noise, 

a  pair  of  neighboring  edge  pixels  parallel  to  the  boundary  both  fail  to 
indicate  edges  as  they  should.  If  there  were  no  correlation  between  neigh¬ 
bors,  P^p  would  simply  equal  Pj^^,  the  square  of  the  probability  of  non¬ 
detection  of  a  single  edge  pixel.  However,  as  we  shall  see,  there  is  a 
positive  correlation  between  neighbors  so  that 


P 

2P  *  N 


2 


(D-1) 


We  specialize  to  neighboring  pixels  along  the  right  boundary  and  cal¬ 
culate  the  (x,  y)  and  (x,  y+1)  correlations  of  the  Sobel  operator  components. 
Rewriting  Eqs.  (7-4)  and  ('’-5)  in  the  main  text  for  the  (x,y  +  l)  pixel,  we 
obtain  tlie  Sobel  components: 


Sp(x,  y+1)  =  I(x-l,y)  +  2I(x,y)  +  l(x+l,y)  -  I(x-l,y+2)  -  21(x,y  +  2) 


-  I(x+1,  y+2) 


(D-2) 


and 


yl  +  2I(x-l,  y+1)  +  l(x-l,  y+2)  -  I(x  +  1,  y) 
-  2Hx+l,y+l)  -  I(x+l,y+2). 


{D-3) 


D-1 


find 


Taking  expectations  of  products  of  adjacent  operator  outputs,  we 


^  (sv(^.  y)  Sy(x, 

y+l)| 

1=  (4AT)^  +  8 

(D-4) 

ff:  (s^j(x,y)  S^(x,y+l)j 

1  =  0 

(D-5) 

y+1)) 

1  =  0 

(D-6) 

<^‘(Sy(x,y)  S^j(x, 

y+l)j 

=  0 

(D-7) 

Introducing  a  more  condensed  notation 

S  ^  S(x,y) 

9 

W  :i  S(x,y+1) 

(D-8) 

U  =  Sj  j(x,  y) 

9 

^  Sy(x,y+1) 

(D-9) 

V  =  Sy(x,  y) 

9 

^  S^(x,  y+1) 

(D-IO) 

have 

s  =  u|  + 

1 V  ! 

,  w  =  1 U  J  +  1 V  J 

(D-11) 

The  pair  probability  we  ' 

wisli 

to  calculate  l^^p  is  given  by 

rT 

/  dS  /  d\V  P(S,  W) 

(D-12) 

O 

•'o 

m 


where  P(S,  W)  is  the  proper  bivariate  density  distribution  function 


P,s,w,  =  fav  fav, 

•'-CO  CO  m 


6  (S  -  U  -  V  )  6(W  -  )  (D-13) 


•  P(U,  U^,  V,  V^) 


Since  U,  U^,  V  and  are  linear  combinations  of  random  Gaussian 
variables.  P(U,  U  ,  V,  V  )  is  a  quadra -variate  Gaussian  density  function. 
In  view  of  (D-4)  -  (D-7), 

-(u^+u+^l/a;^  j 

P(U.  U^,  V.  V^)  .  ^  e  .  - ^ - - 

^”2 


exp 


(V  -  4AT)^  -  2p(V  -  4AT)  (V  -  4AT)  +  (V^  -  4AT)^ 


2?^1  - 


(D-14) 


where 


2  ,  ,  2 
cr  1 2  IT 


lD-15) 


and  tbc  correlation  coefficient 


P 


(D-16) 


U-  5 


Substituting  (D-14)  into  (D-13),  and  (D-13)  in  (D-IZ)  leads  to  a  six¬ 
fold  integration.  Although  it  is  possible  to  perform  most  of  the  integrations, 
we  have  not  succeeded  in  performing  them  all  without  making  unsatisfactory 
app  roximations . 

Consequently,  we  will  approximate  D-13  with  a  bivariate  Gaussian 
distribution: 


P(S,  W  )  = 


,  2 

1  2 

2r 

/i  -  Ps 

exjj 


(S  -  S)^  -  2  Pg(S  -  S)(W  -  S)  +  (W  -  S)^ 


(D-17) 


Note,  that  tliis  reduces  to  the  ordinary  Gaussian  function  used  previously, 

—  2 

The  average  intensity  S  and  variance  have  been  calculated 
previously  in  Eqs.  (C-5)  and  (C-11),  respectively.  The  correlation  coeffi¬ 
cient  p  is  calculated  from 
^  s 


Pg  (7,^  =  <f  (S\V)  - 


(D-18) 


/  (S\V)  =  2 


(|u|  •  |u_^  |)  +  <^  (|u|  •  |v^  |)  +  f^(|v|  •  |u_^  |)  +  |v|  •  |v^|) 


(D-19) 


Since  we  are  assuming  a  sufficiently  higli  SNR  as  in  Appendix  C, 


hav( 


^  ( 

|u|  . 

1  V 

1  + 

[)  4AT 

U|)  =  . 

{ 

Ivl  • 

iu 

j )  ^  4AT 

/  ( i 

LI  I  >  =  4 

1  1- 

i 

-t  I 

^  ( 

Ivl- 

I  V 

I  H 

)  'MV 

=  (4AT) 

■VI 


(D-20) 

(D-21) 

(D-Z2) 


D-4 


Since  U  and  are  independent  variables 

<'(|U|  •  |U^|)  =  rf  (|U|)  •  *f  (|U^|)  =  (D-Zi) 

Solving;  Eqs.  (D-14)  -  (D-19)  for  p^with  the  aid  of  (C-5)  and  (C-11), 
one  obtains 


P 


S 


0.  73p 


(D-Z4; 


Since  by  (D-16), 


p=  2/3. 
p  =  0.49 


(D-Z5) 


Tl\e  formal  approximate  solution  to  the  problem  is  now  given  by 
(D-1  Z)  with  tlie  lower  limits  extended  to  -oo.  consistent  wdth  our  assumption 
that  use  of  the  bivariate  Gaussian  distribution  (D-17)  is  justified, 


ZP 


(D-26) 


Introducing  new'  integration  variables  x  and  y  defined  by 


X 


(D-27) 


S  -  W 


(D-28) 


1)  -  4 


obtjui.') 


P 


ZP 


L 


T 


1 


l  .o  i)i‘(.)bability  function  defined  by 


I  (h,k,p) 


dy  g(x,y,p) 


(D-29) 


(D-30) 


■  ui'i'c  y,p)  is  the  bivciriate  normal  probability  function  defined  by 


h(x,  y, p) 


1 


exp 


r 


1 

2 


/x^  -  ^p-^y  + 
\  ,  2 


(D-31) 


The  functioiis  L(h,k,p)  and  e(x,y,  (')  are  defined  in  Section  26,3  of 
Abramo'A'itz  and  Stegun  NBS  Handbook  of  Matliernatical  functions.  Using 
another  fornuila  given  there,  Eq.  (D-29)  can  be  put  in  the  form 


P 


2P 


(D-32) 


P,|,  is  then  obtaijied  by  gi-apliical  interpolation  from  iso-valued  cur\’es  given 
in  t!iis  lu'ference  for  I.(h,o,i))  as  function  of  li  and  f. 


D-u 
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APPENDIX  E.  PROBABILITY  OF  PAIR  GAP 
PERPENDICULAR  TO  BOUNDARY 


In  this  appendix,  we  calculate  the  probability  P^g  that,  due  to  noise, 
a  pair  of  neighboring  edge  pixels  perpendicular  to  the  boundary  both  fail  to 
indicate  edges  as  they  should. 

The  analysis  proceeds  as  in  Appendix  D  except  that  we  are  now  deeding 
with  correlation  along  the  right  boundary  between  pixels  (x, y)  and  (x+l,y). 
Rewriting  the  Sobel  equations  for  the  (x+l,y)  pixel,  one  obtains 

=  I(x,y-1)  +  2I(x+l,y-l)  +  I(x+2,y-l)  -  I(x,y+1) 

-  2I(x+l,y+l)  +  I(x+2,  y+1)  (£-1) 

S^(x+l,y)  =  I(x,y-1)  +  2I(x,y)  +I(x,y+1)  -  I(x+2,y-l) 

-  2I(x+2,y)  -  I(x+2,y+l)  (E-2) 

Taking  expectations  of  adjacent  operator  outputs,  we  find 

^(Sy(x,  y)  Sy(x+l,y))  =  (4AT)^  (E-3) 

<^|Sy.(x,y)  Sjj(x+l,y)j  =  0  (E-4) 

^(Sj^(x,  y)  Sv(x+l,y))  =0  (E-5a) 

^(Sj^(x,  y)  SH(x+l,y))  =  8cr^  =  pa^  (E-5b) 

E-1 


where  as  in  Appendix  D. 


2 

P  =  3  . 

(E-6) 

We  again  assume  the  bivariate  Gaussian  distribution  (D-17)  with  S 
and  W  given  by  Eqs.  (D-11)  with  the  new  definitions 

S  =  S(x,y)  ,  W  s  S(x+l,y) 

(Er7) 

U  =  S^{x,y)  ,  =  Sj.j(x+l,y) 

(E-8) 

V  =  Sy{x,y)  ,  H  Sy{x+l,y) 

(E-9) 

Corresponding  to  (D-18),  the  correlation  coefficient 

is  calculated  from 

Pg  <^2  =  • 

(E-10) 

Eq.  (D-19)  is  still  valid,  but  the  terms  have  to  be  recalculated  since 
the  definitions  and  are  different  than  they  were.  Again  assuming  a 

sufficiently  high  SNR,  we  have 

•  |V^  |)  ^  <?(V  V_^)  =  (4AT)^ 

(E-11) 

by  Eq.  (E-3) 

<^^(|vl  •  |U^|)  =  4AT  <g'(|u^|)  =  4AT^£ 

(E-12) 

^  (|U|  •  |V^|)  ^4AT  <g'(|u|)  =  4AT^a 

(E-13) 

The  calculation  of  <5“(|u|  •  is  more  difficult  since  U  and 

are  correlated,  and  both  have  zero  mean  so  that  the  high  SNR  approximation 
is  no  help  here. 


We  have 


/«  -00 

dU  /  dU^  •  |u|  •  |u^|  •  P(U,U^, 

00  *' •-  CO 


p) 


where  P(U,U^,  P)  is  the  bivariate  Gaussian  distribution  defined  by 

1 


P(x,y,  P)  = 


_  p' 


exp 


2  ^  2 
X  -  2pxv  +  V 

->/!  2,  2 

2(1  -p  )  £ 


(E-14) 


(E-15) 


Integrating  over  each  quadrant  separately,  changing  some  integration 
variables  to  their  negatives,  and  combining  terms,  one  finds 


o5“(|u|  •  |u^|)  =  2  y  dx  /  dy  xyP  (x,  y. 


p) 


r  r 

+  2  I  dx  /  dy  xyP(x,y, -p). 


(E.16) 


We  integrate  the  first  term  first  by  converting  to  polar  coordinates. 


x  =  rcosB,  y  =  r  sin  6  . 


(E-17) 


j\.  f 


dy  xyP(x,  y,  p)  = 


-  P 


f 


sin  2ede  •  F(e)  (E-18) 


where 


■/ 


F(0 )  =  J  r"'  dr  exp 
o 


r  (1  -  p  sinZ6) 
2(1  -P^)£^ 


(E-19) 


(g  -p^)rM 

\  1  -  p  sin  20/ 


E-3 


Therefore 


/. 


sin  zede  •  F(e)  =  2(1  -  ^ 


L 


sin  26d6 


(1  -  p  sin  20) 


2(1  - 


~  f 

0P  / 

n 


(E-20) 


de 


1  -  p  sin  20 


From  integral  tables,  one  finds 


/, 


d0 


1  -  p  sin  20' 


Vl  -  P^' 


^  .  -1 
+  tan 


V7772 


(E-21) 


The  derivative  of  (E-^,!)  with  respect  to  p  is 


n  2. -3/2 

p (1  -  p  ) 


f  +  tan-‘ - E 


V777 


(E-22) 


1  -p 


Substituting  (E-22)  into  (E-20),  and  (E-20)  into  (E-18)  one  finds 


/.dK  I 


ydy  P(x,y,  p)  =  = 


-1 


+  V1T72 


{E-23) 

Substituting  (E-23)  and  a  similar  quantity  with  p  replaced  by  -p  one 


obtains 


(lul  -  lull  = 

0  tan~V  — — 2 \  +  Vl  -  p  ^ 

\l  1  1+1/  TT 

\  \/t  2  j 

V  \1  -  p  / 

(E-24) 


E-4 


I 


Combining  Eqs.  (E-10),  (D-19),  (E-11),  (E-12),  (E-13),  (E-24), 
(C-5),  and  (C-11),  and  solving  for  one  obtains 


o  =  ptan"Mp/(l  -p¥^^)  +  ,£._2  j 

Numerical  calculation  verifies  that  p  is  a  monotonicailly  increasing 

8 

function  of  p  as,  of  course,  it  must  be.  For  the  case  of  interest,  p  =  2/3, 

p^^O.ll.  (E-26) 

s 

Note  that  for  edge  operators  other  than  the  Sobel,  Eq.  (E-25)  for  p^ 
and  Eq.  (D-24)  for  pg  will  still  be  valid  with  different  vadues  of  p. 

The  formal  solution  for  P2S*  probability  of  a  pair  gap  perpendicu¬ 

lar  to  the  boundary,  is  the  same  as  that  given  in  Appendix  D  for  P2P’ 
However,  for  the  value  of  Pg  given  in  (E-26),  the  correlation  is  suffi¬ 
ciently  small  that  neglecting  the  correlation  altogether  introduce  little 
error.  Hence 


P  ~  P 
^2S  =  N  » 


(E-27) 


where 


is  the  probability  of  not  detecting  any  specified  edge  pixel  (excluding 
corners,  of  course). 


(E-28) 


APPENDIX  F.  NEXT-NEAREST-NEIGHBOR  CORRELATION 


In  order  to  gain  some  idea  of  how  good  our  neglect  of  next-nearest- 
neighbor  correlation  is  in  our  theory  of  gaps  higher  than  second-order,  we 
calculate  the  correlation  coefficient  of  the  Sobel  operator  along  the  right 
boundary  between  pixels  (x,y)  and  (x,y+2). 

Rewriting  the  Sobel  equations  for  the  (x,  y+2)  pixel,  one  obtains 

Sj^(x,y+2)  =  l(x-l,y+l)  +  2I(x,y+l)  +  I(x+l,y+l)  -  I(x-l,y+3) 

-  2I(x,y+3)  +  I(x+l,y+3)  (F-1) 

Sy(x,y+2)  =  I(x-l,y+l)  +  2I(x-l,y+2)  +l(x-l,y+3)  -  I(x+l,y+l) 

-  2I(x+l,y+2)  -  I{x+l,y+3)  (F-2) 

Taking  expectations  of  operator  outputs  separated  by  two  pixels, 

we  find 


<^(sy(x,y)  Sy(x,y+2)j 

=  (4AT)^  +  2 
=  (4AT)^  +  py^^ 

(F-3) 

«?|Sy(x,  y)  Sj^(x,  y+2)| 

=  0 

(F-4) 

<g’^Sjj(x,y)  Sy(x,y+2)| 

=  0 

(F-5) 

<^(sH(*.y)  Sjj(x.y+2)) 

4  2  „  2 
=  -6«r  =  ppj?  , 

(F.6) 

F-l 

where 


Py  =  1/6  (F-7) 

and 

=  -1/2.  (F-8) 

V\  e  again  assiune  the  bivariate  Gaussian  distribution  (D17)  with  S 
and  W  given  by  Eqs.  (D-H)  with  replaced  by  the  next-nearest- 

neighbor  correlation  coefficient  which  we  wish  to  calculate,  and  (D-8), 

(D-9),  and  (D-10)  replaced  by 


S  =  S(x,y)  ,  W  =S(x,y+2) 

(F-9) 

U  E  Sj^(x,y),  =S^(x,y+2) 

(F-10) 

V  E  Sy(x,y),  E  S^(x,y+2) 

(F-11) 

Corresponding  to  (D-18),  the  correlation  coefficient  is  calcu¬ 

lated  from 

PnNN  '2^  =  «®W)  -  S-2 

(F-12) 

Eq.  (D-19)  is  still  valid,  but  the  terms  have  to  be  recaJcuIated 
the  definitions  of  and  have  changed.  Again  assuming  a 

high  SNR,  we  have 

,  since 
sufficiently 

<?(lv|  .  |v^|j  ^  <?(VV^)  =  (4AT)2  + 

(F-13) 

by~Eq.  (F-3). 


F-2 


I 


'^(KI  •  •^(l“+l)= 

#(|U|  .  |VJ)=.4AT  ,f(|u|)  = 

The  calculation  of  «f(|u|  *  [u^|)  is  given  by  Eq.  (E-24)  with  p 
replaced  by  Combining  all  the  necessary  equations  as  previously, 

obtains 


^NNN 


1T 


W'  -  Oh" 


-  1 


V 


Applying  the  values  for  Py  and  pj^  given  in  Eqs.  (F-7)  eind  (F- 
one  obtains 


This  amount  of  correlation  is  still  sufficiently  small  that  it  is  a 
satisfactory  approximation  to  neglect  it. 


(F-14) 

(F-15) 

one 

(F-16) 

8), 

(F-17) 


F-3 


APPENDIX  G.  DERIVATION  OF  FISHER  INFORMATION  MATRIX 

FOR  GENERAL  POLYGON 


In  this  appendix  the  basic  derivations  necessary  for  the  results  of 
Section  4  are  presented. 

Since  S  is  constant  except  near  the  boundary,  the  Fisher  Information 
Matrix  may  be  written 


T  _  _2_  [ f  9S(x,y,i)  9S(x,  y,i) 

11^.. 


dxdy 


(G-1) 


where  J  Jn  •••  denotes  integration  over  a  narrow  strip  of  width  6 
^  n 

surrounding  the  nth  side. 

In  order  to  more  conveniently  perform  the  integration  and  express 
the  function  S(x,  y,  f)  near  the  nth  side,  we  translate  and  rotate  to  new 
coordinates  U  and  V  (see  figure)  defined  by 


X  =  X  ,  +  U  sin  6  +  V  cos  6  (G-2a) 

n-1  n  n 

y  =  y  1  -  U  cos  0  +  V  sin  0  (G-2b) 

’  'n-1  n  n  '  ' 

u  =  (x  -  sin  -  (y  -  Yn-l^  ®n  (G-3a) 

V  =  (x  -  cos  8,^  +  (y  -  ®n  (G-3b) 


G-1 


2 


The  condition  that  S(x,  y,  f)  decreases  linearly  to  zero  in  the  boundary 
regions  is  expressed  by  the  equation 


S(x,  y,  f)  -  AT 


l.U<-| 

0.  u>| 


So  in  terms  of  x  and  y  (for  points  in  R^) 


(G-4) 


S(x,  y,  f)  =  AT 


1  1  ■ 
2  ■  6 


(x  -  sin  -  (y  -  ,)  cos  6 


n-r 


(G-5) 


Since 


f  =  at  is  _  i_y 

1  ^  '  afj  2  6  ' 


while  for  i  >  1 
ax 


as 

af, 


AT 


ay 


1  '' 7 «  1  rtbn  .  V  •  ooxi 

-  -a-f^"*"  ®n  +  -aV  ®n  +  -  "n-l>  ®n  "af-  ^  ’  ^n-l^  ®n  afj 


aen 

afj 


(G-6a) 


aen 


AT 

6 


ax 


n-l 


aii 


sin  0  + 


ay 


n-l 


af; 


cos  0  +  V 


aen 

af; 


(G-6b) 


where  use  was  made  of  Eqs.  (G-2). 


It  is  straightforward  to  verify  that 


The  integral  in  (G-1) 


,9  ^ 


II 

R 

n 


9S(x,  y,  f)  aS(x,  y,  ^ 


af 


af 


dxdy 


(G-14) 


can  now  be  evaluated  in  terms  of  the  U,  V  variables: 


.9  = 


/;  -  / 


n 


dv 


bs  as 

afi  afj 


(G-15) 


Substituting  (G-13)  into  (G-15),  the  integrations  are  easy  to  carry  out 
so  to  obtain 


Sx  9y 

®n  3^  ®n  Bi~ 


3x 


“aT"  ^  ®n  9fj 


f 


L^(AT) 
"  36 


-  sin  e 


1 

+  cos  e 

n  9f^  n  9fj 


9^  1  9y^  1 

sin  e_  — +  cos  6. 


n  9f. 


"n  afj 


1 

^  2 


9x 


9yr 


®n  ^  ^n 


i  J 


9  X 


sin  e 


n-1 


9y 


n  9f . 
3 


+  cos  e 


n-1 


n  9f. 
3 


4 


a>^n-l  ^^n-l 

-  sin  6  — 57 —  +  cos  e„  sj — 

n  9f.  n  91. 

1  1 


9x 


3y. 


-  sin  6  +  cos  0  -g-r- 

n  9;.  n  di. 

3  3 


(G-16) 


Introducing  the  midpoint  variables  x^,  y^. 


X  = 

n 


X  1  + 
n-1  n 


(G-17) 


^n  = 


^n-l  ^n 


(G-18) 


G-4 


the  bracketed  ejqiression  in  (G-16)  simplifies  to 


ax 


n 


ay 


-  sin  G„  ^7^  +  cos  e 
n  di. 

1 


n 


n  af. 


ax 


-  ®n  -afT  +  ®n  af. 

J  J 


aAx 


-  sin  6 


n 


aAy 


n  af. 


+  cos  0 


n 


n  af. 

1 


aAx 


-  sin  6 


n  af 


+  cos  6 


aAy^ 


J 


n  af 


J 


(G-19) 

Combining  (G-19),  (G-16),  (G-15),  (G-14),  (G-7b)  and  (G-1)  we  obtain 
N 


J  =  2(AT)^  Y' 

ij  ''  No 6  Z— I 


n=l 


n 


dx 

-  sin  6  -^7^  +  cos  6 
n  df. 


n  af. 

1 


dx 


dx 


^  n  .  ^  n 

-  sin  0  -^7-  +  cos  0 


n  af. 


n  afj 


L  2  ) 

n  a0n  a6n  ( 

12  at.  atj  j 


(G-20) 


The  quantities  J^j  are  similarly  obtained  using  (G-6a)  and  (G-13)  in  (G-15), 
to  obtain 


= 


^  V  T  I  ^  ^^n-i  «  aen 

No  2-(  ^n  \~  af.  ®n  af.  ®n)'^  2  af. 

n=  1  '  ^  J 


N 

No  /  J  n 
n=l 


ax  a  y 

-  sin  0  -57—  +  cos  0 

n  af.  n  af. 


(G-21a) 

(G-21b) 


N 


AT 

No 


n=l  L 


9y, 


a  X 


»  n  .  n 

"^^n  af.  ■  ‘^^'n  af. 

J  J 


(G-21c) 


G-5 


The  quantity  is  obtained  directly  from  the  original  Fisher 
integral  since  9S/ 9AT  =  1  over  most  of  the  interior.  Thus,  to  a  good 
approximation,  =  2/N^-(area  enclosed  by  boundary).  (G-22) 


G-6 


APPENDIX  H.  EFFECT  OF  CLUTTER  ON  TARGET  SEGMENTATION 

PROBABILITY  OF  ANOMALY 


In  this  appendix,  the  bounds  on  probability  of  anomaly,  i.e.,  the 
probability  that  a  target  segmentor  will  segment  a  candidate  object  nowhere 
near  an  actual  target,  will  be  derived  for  the  case  where  there  is  a  single 
clutter  object  in  the  scene.  This  extends  results  previously  obtained  in  (1) 
for  rectangular  targets  imbedded  in  white  noise.  The  performance  bounds 
on  probability  of  anomaly  will  be  derived  for  three  cases.  The  first  case 
has  a  target  and  a  clutter  object  in  the  field  of  view,  and  the  segmentor  has 
been  optimized  as  if  it  were  a  known  shape  clutter  object  in  the  scene,  but 
the  size  and  location  is  unknown.  In  the  second  case,  the  bounds  will  be 
computed  where  there  is  a  target  aiid  clutter  object  in  the  field  of  view,  but 
the  segmentor  is  not  optimized  for  the  clutter,  i.e.,  no  clutter  is  assumed 
when  the  optimum  segmentor  is  derived.  The  third  case  is  the  reverse  of 
the  second  case  where  there  is  no  clutter,  but  the  target  segmentor  is  opti¬ 
mized  as  if  there  were. 

The  first  model  is  governed  by  the  equation 

I(x,y)  =  T^  S(x,  y.  Fg)  +  T^  (l  -  S{x,  y,  F^  ))  C(x,  y,  F^) 

+  n(x,y)  (H-1) 

where 

•  I(x,  y)  is  the  observed  intensity  distribution 

•  T,j,  is  the  target's  d.c.  intensity  offset 

•  T^  is  the  clutter  object's  d.c.  intensity  offset 

•  Fg  is  a  feature  vector  which  specifies  the  size  and  location  of  the 
target 


H-1 


•  Fj,  is  a  feature  vector  which  specifies  the  size  and  location  of  the 
clutter  object 

j  1  if  X,  y  is  on  the  target 

•  S(x,  y,  Fg)  =  ' 

I  0  otherwise 

I  1  if  X,  y  is  on  the  clutter  object 

•  C(x,  y,  F^)  = 

I  0  otherwise 

•  n(x,  y)  is  a  zero-mean  Gaussian  process  with  autocovariance 
N  /26(t  ,  T  ). 

O'  X*  y 

The  shape  of  the  target  and  the  clutter  object  are  assumed  known, 
but  the  size,  location,  and  d.c.,  intensity  offset  is  assumed  unknown. 

With  these  assumptions  the  optimal  target  segmentation  will  be 
obtained  by  the  method  of  maximum  likelihood.  The  log -likelihood  of  the 
observed  intensity  distribution  is 


a(t,, 


^  ^s(x,  y,  y  (l(x,y)  .-f)  dxdy 
*  //  ('  -  S**'  V.  i,))  C(x,  y,  Q 


e  X 

c 


(l  (x,y) 


dxdy 


(H-2) 


The  values  of  T  ,  T  ,  f  and  f  that  maximize  this  function  are  the 
s’  c  — s  "c 

estimated  parameters.  Noting  that  the  maximum  likelihood  estimates  of 
T„  and  T„  are 


T 


S 


II 


S(x,  y,  f^) 


l(x,  y)  dxdy 


S(x,  y,  fg)  dxdy 


(H-3a) 


H-Z 


and 


jj (l  -  S(x,  y,  £g)|  C(x,  y,  f^)  l(x,y)  dxdy 
jj -  S(x,  y,  fg)|  C(x,  y,  f^)  dxdy 


fche  result  is 


A  =  A 


('^S*  ^C*  ^s'  ic)  “  1 


where 


a2  -  _2- 

1  ~  N 


(Ih 


y)  S{x,  y,  f  )  dxdy 


JJs(k,  y, 


)  dxdy 


and 


.2  2 


t 

(//<■  -  S(x,  y,  f^)j  C{x,  y,  f^)  I(x,  y)  dxdyj 


(H-3b) 


(H-4a) 


(H-4b) 


(H.4c) 


JJ{^  '  y.  £g))  C{x,  y,  f^)  dxdy 

In  order  to  compute  probability  of  anomaly  it  should  be  noted  that  A  j 


and  A^  have  the  following  properties: 


1.  Aj^  and  A^  are  statistically  independent 


2.  E(A,) 


) 


H-3 


3.  EiA^) 


(1  -  S)  CS^  +  T 


(1  -  S)  C  (1 


S^) 


S)C 


4.  Var  Aj  =  1,  Var  A^  =  1 

5.  Aj  and  are  Gaussian  random  variables. 

(The  notation  and  is  used  abov'e  to  denote  the  functions  S(x,  y,  F  ) 
and  C(x,  y,  F^)  when  Fs  and  Fc  are  equal  to  the  target's  and  clutter's  true 
featr.re  vector  Nalues,  and  the  notation  S  and  C  is  used  to  denote  the  func¬ 
tions  S(x,  y,  Fs)  and  C(x,  y,  Fc)  with  Fs  eUid  Fc  equal  to  the  hypothesized 
values  of  target's  and  clutter's  feature  vectors.  This  notation  will  be  used 
throughout  this  appendix.  All  integrations  are  performed  over  the  image 
plane.  ) 

Using  the  five  properties  of  Aj  and  A2  given  above  and  a  well  known 

relationship  between  the  sum  of  squares  of  Gaussian  random  variables, 

2  2 

the  probability  distribution  of  A=  A^  +  A^  can  be  determined.  This  rela¬ 
tionship  will  be  stated  in  theorem  form.  Theorem:  Let  x^,  i  =  I,  Z,  3,  ..., 
N  he  independent  Gaussian  random  variables  with  zero  mean  and  variance 
equal  to  1,  Then 


irr  1 


2 
X  . 
1 


is  a  random  variable  having  the  non-central  chi-square  distribution,  X  (N,  a), 
with  N  degrees  of  freedom  and  nun-centralit>  parameter 


The  cumulative  distribution  function  (c.d.f. )  of  the  non -central  chi-square 
with  N  degrees  of  freedom  is  given  by 


P(Q  s  t)  =  ^  ^ 

j=0 


whe  re 


1  N 


=  [2«/^-(N/2)]  r  y'"' 

J  o 


e-y/2  dy 


is  the  c.d.f.  of  the  central  chi-square  random  variable. 

Using  the  five  properties  of  A  j,  and  and  the  theorem  that  has 
been  just  stated,  it  is  seen  that 

2  2 

A=  Aj  +  A2 


has  a  non -central  chi-square  distribution  with  two  degrees  of  freedom  and 
non-centrality  parameter 


H-5 


For  S  =  S,p  and  C  =  C^,  i.e.,  the  hypothesized  target  and  clutter  feature 
vectors  values  are  equal  to  the  true  target  and  clutter  feature  vectors 
values. 


(H-6) 


As  done  in  reference  1,  we  will  define  the  "anonrialy  threshold"  as 


\  = 


C  = 


) 


(H-7) 


and  so 


X  =  E  (x^(2.o))  =2+0 


By  definition,  an  zinomaly  has  occurred  if  A  >  A  for  JJ  SS^  s  0,  that  is, 
an  anomaly  has  occurred  if  A  >  A  for  a  hypothesized  target  region  which  does 
not  intersect  the  true  target  region.  To  compute  the  probability  of  this 
occurrence,  we  return  to  H-5  and  note  that  for  //  SS^  =  0 

'"IF 

(H-8) 


a  = 


No 
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If  it  is  further  assumed  that  the  clutter  is  not  masked  by  the  target  and  that 
this  information  is  available  to  the  target/clutter  detector,  then; 


0. 


Hence, 


Eq.  (H-8)  reduces  to 


a 


N 

o 


(H-9) 


There  are  five  cases  that  the  hypothesized  target  and  clutter  feature 
vectors,  f^  and  will  be  considered  to  be  in  that  may  result  in  an  anomaly. 
Case  (1)  is  defined  such  that  the  hypothesized  target  and  clutter  regions  are 
contained  in  the  region  which  contains  only  noise.  That  is,  for  case  (1) 


jjsc,-  jjs,c-  jj 


cc^  =  0. 


For  any  f^  in  case  (l)  a  =  Oj  =  0.  Hence  the  probability  of  an  anomaly 
for  a  fixed  £^,  f^  in  case  (l)  is  given  by 

P(x^(2,0)  >  k) 

For  all  hypothesized  target  regions  of  area  A(S)  and  hypothesized  clutter 
regions  of  area  A(C)  the  probability  of  anomaly  due  to  region  (1)  can  be 
approximated  by 

(‘  (A(gWA<c)  ) 

where  the  scene  is  assumed  to  have  dimension,  R  x  R. 


H-7 


A  bound  on  probability  of  anomaly  for  all  sizes  of  hypothesized  target 
and  clutter  regions  can  be  obtained  by  applying  the  argument  following 
Theorem  1  in  reference  [l].  The  probability  of  anomaly  due  to  region  (1)  is 
then  seen  to  be  bounded  by 


^  ITT - — - \  -  ^2 

2Ma(S^.^)  +  McA  \  '  / 


(H-16) 


or,  equivalently. 


A(S  .  )  +  A(C  .  ) 
min  mm 


(H-17) 


where  A(S  .  )  and  A(C  .  )  are  the  smallest  hypothesized  target  and 
'  min'  min'  ® 

clutter  areas  allowed. 

Case  (2)  is  defined  such  that  the  hypothesized  target  region  does 
not  intersect  the  true  clutter  region  and  that  the  hypothesized  clutter  region 
is  contained  in  the  true  clutter  region.  For  f  and  f  in  case  (2)  Eq.  (H-9) 
reduces  to 


“  =  ^2  = 


ih 


(H.18) 


where  C  s  C,j,.  We  will  assume  that  if  an  anomaly  occurs  in  region  (2)  it 


will  be  due  to  some  C  £  such  that 


(//”j 

Ih 


IS  maximum. 


H.8 


and  so 


a 


2 


(H-19) 


Since  there  is  only  one  way  to  choose  0  s  C  s  C^,  such  that 


2 


is  maximum,  the  probability  of  anomaly  due  to  region  (2)  is  approximated 
by 


where 


Case  (3)  is  defined  such  that  the  h^^pothesized  clutter  region  is  contained  in 
the  true  target  region  and  the  hypothesize  target  region  is  contained  in  the 
region  which  contains  only  noise.  From  Eq.  (11-9), 


a  = 


H-9 


1 


where  0  <  C  <  S~,  Then  for  fixed  f  ,  f  in  region  (3)  the  probability  of 
anomaly  is  given  by 


1 


(A) 


Assume  that  if  an  anomaly  occurs  in  region  (3)  it  is  due  to  some 
C  <  such  that 


(H-22) 


is  maximum. 

Using  the  notation  to  denote  a  C  such  that  (H-22)  is  maximum, 
the  probability  of  anomaly  due  to  some  f^,  f^  in  region  (3)  is  approximated 
by 


N 


3 


(H-23) 


H-10 


.  ^  that  the  hypothesized  target  region  is  contained  in 

Case  (4)  is  defined  such  that  th  yp^  ^  ^  contained  in  the 


the  true  clutter  region 


and  the  hypothesized  clutter  region 


region  which  contains  only  noise 


.  For  Eq.  (H-9), 


0=0. 


L 


If 


c^s 


<  S  £  c^.  Then  for  any  fixed  f^,  f^  m  case 


(4)  the  probability  of 


where  0 
anomaly  is  given  by 


where 


S  <  C^.  Assuming 

S  <  C,j,  such  that 


that  if  an  anomaly  occurs 


due  to  region  (4)  is  due  to  some 


2 


S 


H-11 


is  maximum,  approximation  of  the  probability  of  anomaly  due  to  region  (4) 
is  given  by 


whe  re 


4  N 


CtS»)  tJ 

IW 


(H-E4) 


N4  = 


UK 

UN 


Case  (5)  is  defined  such  that  the  hypothesized  clutter  region  is  contained 
in  the  true  target  region  and  the  hypothesized  target  region  is  contained  in 
the  true  clutter  region.  For  fixed  f^,  f^  in  case  (5)  the  probability  of 
anomaly  is  given  by 


1  -  F,  (A) 


whe  re 


It  //cSt 

Ih 


(/M 

//= 


H-12 


Using  tlie  results  from  cases  (3)  and  (4),  the  approximate  probability  of 
anomaly  due  to  case  (5)  is  given  by 


maxINj,  N^).  (l  -  (H-25) 

Combining  equations  (H-17),  (H-20),  (H-23),  (H-24)  and  (H-25)  the  bound 
on  probability  of  anomaly  is  seen  to  be 


P(anomaly)<  (ijs".  f+\(c  .  )) 

\  '  min  '  min  /  \  '  / 

,  (l  -  t  N^jl  -  F,^^^(X)j  (H-26) 

where 

"2  j/^T  =  "  •  •  jj^T 

CNR  is  the  "clutter  to  noise  ratio" 


SNR  is  the  "signal  to  noise  ratio" 


N 


3 


a 


4 


N 


4 


Now  consider  the  case  '.'.here  there  is  a  target  and  a  clutter  object  and  the 
segmentator  has  been  optimi/ed  for  only  a  target.  The  livelihood  for  this 
case  is 


A  = 


(// 


I(>^»  y)  S(x,  y)  dxdy 


//s(x. 


S(x,  y)  dxdy 


(H-26b) 
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From  Eq.  (H-26)  the  "ainomaly  threshold"  is  seen  to  be 


A  =  ^  (A|S  =  S„)  = 


There  are  two  ways  that  an  anomaly  can  occur.  The  target  detector  coiild 
lock  onto  noise,  or  the  target  detector  can  lock  on  to  the  clutter  object. 

The  probability  that  the  target  detector  locks  onto  noise  is  obtained  directly 
from  reference  1: 


P(anomaly  due  to  target)  u 


detector  locking  onto  ^ 


A(S  .  ) 
min 


(4k) 


A(S  .  ) 
'  min 


SNR^  •  A(S^)  + 


The  probability  that  the  target  detector  locks  onto  the  clutter  object  can  be 
approximated  by 


fh 

//■• 


P  (  A  >  \  S  =  S,.,)  = 


Ih 

II- 


(l  -  srf'.  (- 


+  «£c.  (  sNr2  +  .  2CNR' 
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Then  probability  of  anomaly  is 


Now  consider  the  case  where  the  segmentor  has  been  optimized  for  a  target 
and  a  clutter  object,  when  there  is  no  clutter.  Probability  of  anomaly  for  this 
case  cein  be  approximated  by  setting  =  0,  =  0  in  Eq,  (H-Z6).  The 

probability  of  anomaly  is  then  seen  to  be 


where 


H-16 
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